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A molecular model developed for resonances observed in medium light heavy-ion colli- 
er. ■ sions is described. At high spins in 2S Si + 28 Si (oblate- oblate system), a stable dinuclear 
configuration is found to be equator-equator touching one. The normal modes around the 
equilibrium are investigated. These modes are expected to be the origin of a large number 
of resonances observed. Furthermore, due to the axially asymmetric shape of the stable con- 
figuration of 2S Si + 28 Si, the system rotates preferentially around the axis with the largest 
moment of inertia, which gives rise to wobbling motion (jf-mixing). Energy spectra for the 
normal modes and for the extended model including the wobbling motion are given. 
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§1. Introduction 



^ I Intermediate resonances observed in heavy-ion scattering have offered intriguing 

subjects in nuclear physics. High-spin resonances well above the Coulomb barrier in 
the 24 Mg + 24 Mg and 28 Si + 28 Si systems exhibit very narrow widths, which suggest 
rather long lived compound nuclear states 

Betts et al. firstly observed a series of resonance-like enhancements at 9 cm = 90° 
in elastic scattering of 28 Si + 28 Si, in the energy range from Ey^ = lOlMeV to 
128MeV with broad bumps of about 2MeV width. They gave spin assignments of 
J = 34 ~ 42 by the Legendre-fits to the elastic angular distributions for each bump, 
■ which correspond to the grazing partial waves.®'® 

| They further closely investigated angle-averaged excitation functions for the elas- 

tic and inelastic scattering in the energy region corresponding to J = 36 — 40, and 
found, in each bump, several sharp peaks correlating among the elastic and inelastic 
channels.®'® The total widths of those resonances are about 150keV, and the inelas- 
tic decay strengths are enhanced and stronger than the elastic one, which suggests 
that they are special eigenstates of the compound system. Similar sharp resonance 
peaks are observed by Zurmiihle et al. in the 24 Mg + 24 Mg system.® The level 
densities observed in those systems are over one per MeV, which suggests activation 
of internal degrees of freedom, in addition to the radial motion. The decay widths 
of the elastic and inelastic channels up to high spin members of the 24 Mg or 28 Si 
ground rotational band exhaust about 30% of the total widths, whereas those into 
a-transfer channels are much smaller.® ™ These enhancements of symmetric-mass 
decays strongly suggest dinuclear molecular configurations for the resonance states. 
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It is also noted that the widths of the elastic channel are rather small, for example, 
a few keV, being quite different from high spin resonances in lighter systems such 
as 12 C + 12 C and 16 O + 16 O, which are well explained by the Band Crossing Model 
(BCM)P^ i.e., by couplings between the relative motion of the incident ions and the 
low- lying collective excitations of the ions. 

From viewpoints of nuclear structure studies, one immediately thinks of sec- 
ondary minima in fission of heavy nuclei, or of superdeformations which have been 
intensively studied in medium mass nuclei.^ Actually Bengtsson et al. made 
Nilsson-Strutinsky calculations for shape isomers of 56 Ni and obtained an energy 
minimum at large deformation, which appears to correspond to a dinuclear con- 
figuration.^ Recently, a couple of microscopic calculations have been performed 
with expectation of existing shape-isomer bandsP^'^ For the ^Mg + ^Mg system 
(A = 24 or 26) with very high spins, theoretical works were made to obtain stable 
configurations.^'^ All those models, however, are not able to reproduce the level 
density of the sharp resonances as well as the decay properties observed. 

Taking into account the difference from resonances in lighter systems and the 
level density of the sharp resonances observed, we have proposed a new dinucleus- 
molecular model for the high spin resonances in the 24 Mg + 24 Mg and 28 Si + 28 Si 
systems P^'^'^ in which two incident ions are supposed to form a united composite 
system. It rotates as a whole in space with the internal degrees of freedom originating 
from interaction of the deformed constituent ions. This is in contrast with the 
viewpoint of BCM. 

Actually, we have already applied the model to the 24 Mg + 24 Mg system to 
obtain a stable dinucleus configuration, using the folding potential. Normal modes 
of motion around the stable minimum were solved with harmonic approximations, 
and several characteristic modes were obtained, such as butterfly one, etc., which 
are expected to be responsible to the observed sharp resonances. Decay properties 
of those resonance states were analyzed and strong enhancements to the mutual 
excitation channels are obtained in agreement with experiments for the 24 Mg + 24 Mg 
system P3 

The same model has been applied to 28 Si + 28 SiP^ As is expected from the 
experience on 

24 Mg + 24 Mg] there 

are several intrinsic modes with excitation energy 
of a few MeV to several MeV. They, thus, are expected to correspond to the sharp 
resonance peaks within each bump of the grazing J. Therefore, the present model 
appears promising also for the sharp high spin resonances observed in 28 Si + 28 Si. 

Recently, a new development has been obtained, giving attention to the remark- 
able difference between the 24 Mg + 24 Mg and 28 Si + 28 Si systems. In the former, 
the stable configuration is pole-to-pole one due to the prolate deformation of 24 Mg, 
while in the latter, it is the equator-to-equator configuration due to the oblate de- 
formation of 28 Si. Therefore, the former composite system is axially symmetric in 
the equilibrium, while the latter is triaxial. Then, in the latter, strong fT-mixing is 
kinematically induced, and results in a wobbling motionPD'GZ) Hence we have ex- 
tended our molecular model so as to include couplings between states with different 
-fC-quantum numbers (projection of the total angular momentum on the molecular 
z'-axis). As a result, we have obtained new low- lying states due to a triaxial shape 
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of the equilibrium configuration. In practice, we do not treat Coriolis terms in the 
hamiltonian explicitly, but we diagonalize the hamiltonian of the asymmetric rotator 
to obtain the rotational spectrum. 

Since the two different kind of models, i.e., the dinuclear molecular model and 
the asymmetric rotator are used to obtain the results, it is necessary to clarify the 
relation between them. We have studied simple examples of dinuclear systems by 
using the molecular model, and have found that the molecular model hamiltonian 
reduces to that of the asymmetric rotator in the sticking limit P3 

The present paper has the twofold aim. One is to describe the molecular normal- 
mode analyses^ as the full paper, and the other is to describe the development 
newly obtained. As for the former, a brief reminder of the molecular model is given 
in §2, where we present the coordinate system and the model hamiltonian in the 
rotating molecular frame. In addition, in Appendix A, we take up simple examples 
of quantization, to compare the kinetic energy expression described by the angular 
momenta in the laboratory frame with that in the molecular frame, and to clarify 
the role of the Coriolis terms. There, the sticking-limit condition of sharing the 
total angular momentum between the orbital motion and the fragment spins is also 
derived. In §3, structures of the 28 Si + 28 Si system are investigated. We begin 
with inspecting the multi-dimensional energy surface and look for the equilibrium 
configuration of the system. In §3.2, harmonic approximation is adopted to solve 
normal modes around the equilibrium. Firstly, an energy spectrum with good K- 
quantum numbers will be given. The symmetries of the system and the practical 
expressions of the wave functions are described in Appendices B and C. 

In order to present the new development, section 4 is devoted to the analyses 
for the dinuclear system with axial asymmetry, which gives rise to wobbling motions 
(i^-mixing) in extremely high spins. After iT-mixing, the instates are recombined 
into new states. The sequence of energy levels obtained by the diagonalization of the 
asymmetric rotator hamiltonian is given. A simple analytic solution is also discussed. 
In §4.2, we take up simple examples of the molecular model hamiltonians and see 
how they reduce to the asymmetric rotator hamiltonians. As a summary, in §5, we 
discuss on the structures of the 

28 Si _ 28 Si 

molecule theoretically explored. 
Those molecular states are expected to be the origin of a large number of res- 
onances observed, and hence theoretical analyses have been made. The results are 
in good agreements with the experimental data,®'^ which will be given in the 
succeeding paper, no. II 

§2. Dinuclear molecular model of the oblate-oblate system 

First, we briefly recapitulate the new molecular model for heavy-ion resonances. 
Definitions and derivations of the expressions are given in detail in Ref. [T5|) , for the 
prolate-prolate system. We have already proposed a new description of interacting 
two oblate-deformed nuclei such as 28 Si + 28 Si.H3'l2Sli j n §2.1, we extend our con- 
sideration to the coordinates of the system including axially-asymmetric deformed 
constituent nuclei, for later descriptions in §4. Subsections 2.2 and 2.3 are devoted 
for the descriptions of the kinetic energy and the nucleus-nucleus potential, respec- 
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tively, some expressions of which are already published in Ref. [T9]) and in a part of 
Ref.[20D. 

The total system is described in terms of rotation of the whole system in space 
and of internal motions of all the other degrees of freedom such as the orientations 
of the deformation axes of two nuclei relative to the rotating molecular axes. We 
anticipate that there exists a stable geometrical configuration, i.e., a minimum in the 
potential energy of the internal degrees of the system. Actually as explained later, 
each typical stable configuration appears by strong attractive nuclear interaction be- 
tween tips of two deformed nuclei, for prolate or oblate deformations, respectively. 
Accordingly, motions of their pole orientations should be treated as vibrational de- 
grees of freedom around the geometrical equilibrium configuration, which is quite 
different from the usual description using "channels" in the weak coupling picture 
with the orbital angular momentum and the spins of the interacting nuclei. 

2.1. Coordinate systems 

The total system to be solved consists of two deformed nuclei interacting with 
each other. We expect the axial symmetry of the constituent nuclei and their con- 
stant deformations, corresponding to the states of the K = ground rotational band. 
We thus start with seven degrees of freedom illustrated in Fig. 1(a), that is, the rela- 
tive vector R = (R, 62, 9\) and the Euler angles of the interacting nuclei (Si, fix) and 
{dt2ih)i where the deformations of the constituent nuclei are taken to be oblate for 
28 Si nuclei.^ When the constituent nuclei contact and interact strongly with each 
other, their deformations in the ground state may change, i.e., additional deforma- 
tions may be induced, such as those associated with the surface 7-vibrations or the 
static asymmetric onesP^ 1 In that case, we have additional degrees of freedom, 71 
and 72, by which the nuclei rotate around their intrinsic z-axes. (Those degrees of 
freedom are not illustrated in Fig. 1, for simplicity.) We define the rotating molecular 
axis z' of the whole system with the direction of the relative vector of two interact- 
ing nuclei, as is shown in Fig. 1(b). In the molecular model, the intrinsic axes of 
each deformed nucleus are referred to the molecular frame as usual. We introduce 
new Euler angles of the interacting nuclei in the molecular frame (ai,/3i,7i) and 
(«2>/52i72) as in Fig. 1(b), which are related to (aj,/3j,7i) by 

Qi (an, fa, 74) = (9x ,0 2 )f2i {pn, A, ji) , i = 1, 2, (2-1) 

where iTs denote Euler rotations with respective angles, with indications of the 
rotations for the each constituent nucleus no. 1 or no. 2 by i. To obtain the configu- 
ration of Fig. 1, it may be more useful to describe in terms of successive rotations as 
i?i(5j, fti, 7j) = Q'i((Xi, j3i, 7i)i?Af(^l) $2)5 where the second rotations Q\ refer to the 
molecular axes, i.e., to be operated on the intrinsic axes which are parallel to the 
rotated molecular axes (x',y',z r ). Correspondingly J2m(#i, #2) appears to be multi- 
plied from the right-hand side in this case. Later, we introduce #3 = (ai + 02) /2 
as the third Euler angle for the rotation of the total system. The axes x' and y' in 
Fig. 1(b) indicate the axes after the rotation Hm{0i,02), while the molecular axes 
x'(6s) and y'[9z) indicate after the whole rotation J?m(#1) #2> ^3)- Note that for the 
constituent nuclei with the axial symmetry, 7, are not necessary, and we put 7$ = 0. 
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Fig. 1. The coordinates of an interacting dinuclear system, (a) shows the relative vector R — 
(R, 62,81) and usual Euler angles (on, Pi) of the i-th nucleus referring to the laboratory frame. 
In (b), the molecular z'-axis and the seven degrees of freedom of the system are displayed, where 
the distance R is not indicated explicitly. The third angle #3 is defined by #3 = (ai + 02) /2 to 
give the whole rotation around the z'-axis. 



Generally we obtain 7j not to be zero due to the transformations between the coordi- 
nate systems, but these 7$ are physically meaningless. They appear in the rotational 
matrices, but they practically disappear in the inertia tensor of the total system; see 
Appendix B of Ref.HB]). 

Large deformations of the constituent nuclei may be induced in the deeply touch- 
ing configurations of the resonances, then their axial symmetry of the deformations 
would be lost and the 7j-degrees of freedom may appear. In the scope of the present 
model, we are able to introduce those degrees of freedom. However without infor- 
mation on the extent of the induced deformations nor on the dynamical properties 
about the deformations in the touching configurations, such efforts would bring no 
fruitful result. Later in §4.1, we consider configurations with such large induced 
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(62,00 



Fig. 2. The dinuclear configuration and the coordinates in the rotating molecular frame for an 
oblate-oblate system. The molecular z'-axis and the seven degrees of freedom of the system 
are displayed, where the ol\- and a2-degrees are to be combined into #3 = (ai + a.2)/2 and the 
degree of twisting a — (ai — 0:2 )/2. The figure is the same as published in Refs. 1 19f> and!20p. 

deformations, and investigate an example of the molecular model with the 7-degree 
of freedom in §4.2, but in this section we restrict ourselves to descriptions without 
induced deformations. This means that we assume the axial symmetry of the con- 
stituent nuclei with their moments of inertia I x = I y and I z = in their principal 
axes, and that we start with the seven degrees of freedom as already mentioned. 
They are illustrated again in Fig. 2; the relative vector (i?, 02,0i) and Euler angles 
of the interacting nuclei in the molecular frame (a%, fix) and (02, ft). The variables 
a± and ct2 are combined into variables #3 = (a± + 012) /2 and a = (ai — 02) /2. Then 
we have 



where Q\ , 62 and #3 are the Euler angles of the rotating molecular frame with the 
other fours being internal variables. 

2.2. Kinetic energy of the dinuclear molecule 

Firstly we obtain an expression of the kinetic energy operator in terms of the 
above coordinates. We start with the classical kinetic energy of the system, which 
can be given in terms of the energies associated with the relative motion (the radial 
motion and the rotational motion of the two-ion centers) and the rotational motions 
of the two constituent nuclei, 



where R denotes the relative distance between the two-ion centers, [i being the 
reduced mass m\m2j '(mi + 7712) of the two nuclei with masses m\ and 777,2, and 
the cm. energy of the total system is omitted. The second term of the r.h.s. of 
Eq. (|2-3p is the rotational energy of the two-ion centers given by the angular velocity 
of the molecular frame u' and the moment of inertia tensor I^R). The diagonal 
components I\\ and I22 of the inertia tensor are (1R 2 , the others being zero, which is 
associated with masses mi and 777.2 at the relative distance R. Then the expression 
of the rotational energy is equal to usual one, \\iR 2 [Q\ + 9f sin 2 62). The vectors u>i 



( qi ) = 2 , 03, «, A, ft) 





(2-3) 
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and ct>2 denote the angular velocities of the rotational motions of the two constituent 
nuclei, t Ui being the transpose of u^. The inertia tensors of the two nuclei I\ and 1 2 
are defined in the coordinate frames of their principal axes. Then, they are diagonal, 
elements of which are determined by the excitation energies of the members of the 
ground rotational bands of the constituent nuclei. 

At this stage, the angular velocities of the constituent nuclei Ui in Eq. (|2-3p 
are still those referred to the laboratory frame, so we have to express them in the 
molecular coordinate system, i.e., in terms of the angular velocity of the molecular 
frame u' and those uj'( referred to the molecular frame. Then we express the total 
kinetic energy as a sum of three parts, the total rotational energy T rot associated 
with a/, the internal kinetic energy 7] n t and the Coriolis coupling term Tq, as follows; 

T = T rot + T int + T c , (2-4) 

T rot = ^Jl s J, (2-5) 

T int = \»R 2 + + ^2 J 2 Cc/ 2 ', (2-6) 

T c = l J {'i?'(ai/3i7i) + *#(a 2 /3 2 72) iW 2 '}, (2-7) 

where R'(ai(3i'ji) denotes the transformation matrix (rotation matrix) which con- 
nects the axes of the molecular frame and the principal axes of each constituent 
nucleus. The total rotational energy T ro t is the rotational energy of the interacting 
constituent nuclei as a whole system, which rotates with the angular velocity u' . 
The inertia tensor is given by 

/ s = I^R) + *i?'(ai/3i7i) Ii#(ai/3i7i) + ^'(^M ^R'i^M, (2-8) 

where the first term of the r.h.s. denotes just the moments of inertia of two-ion 
centers, and the second and third terms are contributions from the constituent nuclei 
individually, though their "rotations" are already taken into account in Eq. (|2-6|) . 
The internal kinetic energy Tj n t is those associated with the orientation degrees of 
freedom of the constituent nuclei in addition to the radial motion between them. 
The last two terms of the r.h.s. of Eq. (|2-6[) have a form of rotational energy, 
but their motions are not necessarily rotational. This is why the quotations are 
put on the word rotations above. Actually the nucleus-nucleus interaction favors 
cohesion of two constituent nuclei, which obstructs rotations of the constituent nuclei. 
Motions in the orientations are, therefore, not necessary to be rotational but would 
be rather confined, such as a sticking of the constituent nuclei and small fluctuations 
thereabout. In the sticking limit, the angular velocities u'- are zero, while they are 
constant in free rotations. We, of course, anticipate intermediate states between the 
sticking limit and the rotation, i.e., fluctuations around the sticking configuration. 
For vibrational motions, for example, we consider fluctuations of the values of u'- 
around zero, average values of them being to be zero. 

After expressing those angular velocities with time derivatives of the correspond- 
ing Euler angles, we obtain a classical kinetic energy expression \ Yl,9ijii^j- And 
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then we quantize it by using the general formula for the curve-linear coordinate 
system, 

(2-9) 

where g and g _1 denote the determinant and the inverse matrix of (gij), respectively. 
As the classical kinetic energy consists of the three parts, i.e., the total rotation, 
the internal motions and their couplings, the quantum mechanical operator for the 
kinetic energy T is also given as a sum of three terms, T = T rot + T[ n t + Tq . Naturally 
the term T rot is associated with the rotational variables {9\,9 2 ,6^), ?l nt with the 
internal variables (R, a, fa, fa) and Tq with both. According to the derivation, T rot 
is expressed by the partial differential operators of 9{. We combine those differential 
operators into angular momentum operators J[ referred to the molecular 
usual, i.e., 

f rot = Y E WW,* ( 2 - 10 ) 
1<«<3 
l<i<3 

where the matrix (i is the submatrix given later, and J-'s are the angular momentum 
operators in terms of the Euler angles of the molecular frame, 



J'l 
J> 2 



cos 6*3 d . d d 

-^ 2 w 1 +sm93 w 2 +cote2COse:i w 3 

sm 9 3 d d . d 

^ 2 w l +cos0 - i w 2 - cotd2Sme3 w 3 



(2-11) 



Here Eq. (|2- 10|> has a form just expected from the classical expression Eq. (|2-5|) . 
but it should be noted that the submatrix fj, is not exactly equal to the inverse of 
the inertia tensor I 8 , due to the Coriolis coupling. The coefficients fn* are given as 
follows, in terms of the internal variables (R, a, fa, fa), 



1 



A*12 
M13 

M23 

^33 



A*22 : 
0, 

1 



fiR 2 ' 



2fiR 2 
1 



cos a (cot fa + cot fa), 



(2-12) 



2fiR 2 



sin a (cot fa 



1 

4L 



1 



+ 



1 



- cot fa), 

1 / 1 

+ 



+ 



2/xi? 2 



I A fiR 2 J sin 2 fa V/ B nR 2 J sin 2 fa\ 2fiR 2 
1 

cos 2a cot fa cot fa , 



where I a and I-q are the diagonal elements of the inertia tensors I\ and I 2 , respec- 
tively. In the definition Eq. (|2-11|) of the total angular momentum operators J- in 
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the body- fixed frame, we write them in terms of the Euler angles 9{, which appear 
at the same time in the coordinates of the relative vector (R, 9, tp) between the two 
constituent nuclei, as 9 = 02 and ip = 9±. Therefore the definition may be misleading 
as not to be the total angular momentum but to be the orbital angular momentum 
L. In Appendix A, we take up simple examples of the quantization both in the 
laboratory frame and in the molecular frame, in order to see the relations between 
the coordinate sets and the definitions for the corresponding angular momentum 
operators. There, the role of the Coriolis coupling term is also clarified. 

The internal kinetic energy operator is associated with the variables (R, a, Pi, fa), 
as already mentioned. As usual, we introduce a volume element dV = dRdadfiidfo 
instead of the original dV = DdRdadfiidfa with D = /j, 3 / 2 R 2 Ia sin PiIb sin/32, which 
means that the wave functions are defined with the additional factor \f~D. Accord- 
ingly we obtain 

Oint + Kidd, 

2 IfidR 2 



'-int 



o mt = -— — 



+ 



+ 



i \ 



fiR 2 Jd(3f 



+ 



(5 



(2-13) 



+ 



i x d 2 
nR 2 ) dp 2 



+ 



2 cos 2a d 2 



fiR 2 dfiidfo 



lf/J_ _L_\ 1_ /J_ _L_\ 1_ 2 \ d 2 

+ 4\ V/a + JJp) ^Wi + V/b + shi 2 /^ ~ ^R 2 J ~cW 



Kdd 



d cos 2a 
da 2fiR 2 

— ^ (cot# 

1 



2fiR 

h 2 



d 

COt Pi COt 02 TT - 

oa 



d n d 

d(3i d/3 2 



sin 2a 



I A + ilR 2 



R 2 )(sin 2 Pi 



+ 1 



)( 



8_ d_ 

da da 



sin 2a 



) 



^ 2 )(sin 2 /3 ? +1 ) 



+ 



2 cos 2a 



fiR 2 



COt Pi COt /?2 



(2-14) 



(2-15) 



where V^dd is the term so-called additional potential due to the new volume element. 

The Coriolis coupling operator Tq consists of coupling operators between the 
variables (#i,#2,#3) and (R, a, Pi, fa), i.e., 



fiR 2 



1 sin a 



( d d \ ±, ( d d \ 

{-m + wJ J ' + cos <dp; + wJ 



d_ 

dfa 



d_ 

dpi 
d 



d_ 

dP2- 



1 
+-sin2a(-cot/3 2 ^- + cot ft — 

2 dPi dp 2 



4 



--(cot Pi 



cot ft 



■(cot ft + cot ft) 



d 9 , 1, d 

— — cos a Ji + Ji cos a—- 

da da 

d . f , f/ . d 
— — sm a Jo + Jo sm a— 

da da 



h 2 



1 



1 



fiR 2 



1 



1 



1 



1 



iiR 2 ) sin 2 ft. 



sin 2 ft 

where the derivative operators of 9i are again rewritten with the angular momentum 
operators ][. 
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For details of some relations and explicit expressions, see Appendices of Ref . fTS]) . 
for example, for the angular velocities in the molecular frame, the classical kinetic 
energy in terms of time derivatives of the Euler angles, their quantization and sym- 
metries of the system. 

In order to make the problem to be tractable, we start with good if-quantum 
numbers firstly, which is expected to be appropriate for the system of small axial 
asymmetry. For the 24 Mg + 24 Mg system (prolate-prolate one), for example, it is 
rather simple to intuitively understand, because stable configurations at high spins 
are dominantly elongated pole-pole ones which keep axial symmetry. However the 
28 Si + 28 Si system (oblate-oblate one) favors equator-equator configurations, which 
do not have the axial symmetry intrinsically. So secondly, the effect of iT-mixing is 
investigated later in §4. 

At this stage, we therefore regroup the kinetic energy operator as follows, 

f =f' + f' c , (2-17) 
f = f; ot + t mt , (2-18) 

where Tq includes the Coriolis coupling Tq and the K- mixing terms in T ro t. Accord- 
ingly the new rotational operator T r ' ot has good i^-quantum numbers. 

Let's restrict our discussion to the rotation and vibration operator T' , together 
with the interaction potential given later. As the kinetic energy operator T' keeps a 
good X-quantum number, eigenstates of the system are of a rotation-vibration type, 

^x- D J MK (6, i ) X K(R,a,f3 u f3 2 ). (2-19) 

Now the problem to be solved is of internal motions, i.e., motions associated with 
the internal variables (R, a, j3\, fa) which couple with each other through the kinetic 
energy operator T' and the interaction potential. For the later use in §3, we define 
the centrifugal potential given by T r ' ot with specified J and K, 



{ J( J + 1) - \k 2 + - cos 2a cot ft cot fa(K 2 - 1)} 



1 1 \ / K z - 1 K- - 1 I 



/ /xi? 2 7 V4sin 2 /3i 4sin 2 /? 2 2 



(2-20) 



where / denotes the moment of inertia of the constituent nuclei, i.e., / = I a = Ib, 
since we are interested in the system of the identical constituent nuclei. In the 
expression of T r ' ot (J, K), we use the eigenvalue K instead of J3. Note that the 
additional potential V^dd in Eq. (|2-15p is moved into T^. ot {J,K) for convenience, and 
similar terms of T r ' ot and T4dd are amalgamated. In the numerical calculations, the 
value of I is estimated from the excitation energy of the 2f state of the 28 Si nucleus. 

2.3. Nucleus-nucleus interaction potential 

For the interaction potential, we want to have an expression that depends on 
geometrical configurations of interacting nuclei, i.e., a potential as a function of the 
Euler angles of the nuclei in addition to the radial distance between them. Proximity 
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potential appears to be one of the most suitable potentials,® but it is rather labo- 
rious to calculate it for various configurations, i.e., one has to find out the shortest 
distance between two curved surfaces of arbitrarily-oriented deformed nuclei and to 
calculate curvatures etc. at the point. Instead, we employ a folding method. Since, 
in the double folding model, nuclear densities corresponding to geometrical molecu- 
lar configurations are directly folded with effective nucleon-nucleon interactions, the 
model easily provides an interaction potential for the present purpose, i.e., as a func- 
tion of the collective variables. As for the nucleon-nucleon interaction, we employ 
one that is called density dependent M3Y(DDM3Y)j^ 

v(E,p,r)=f(E,p)g(E,r), (2-21) 

where f(E,p) gives nucleon-density dependence by 

f(E,p) = C(E)[l + a(E)e-^ E H (2-22) 

p denoting density of nuclear matter in which the interacting nucleons are embedded, 
and g(E, r) describes the original nucleon-nucleon interaction, 



9{E,r) 



e -4r e -2.5r 



7999— 2134- 



4r 2.5r 



+ J{E)5{r). (2-23) 



The first term of g(E, r) is M3Y potential without OPEP and the second term 
represents that from single-nucleon exchange, suggested by Satchler and LoveP^ 1 
E is the bombarding energy per nucleon, which is chosen to be as suitable for the 
resonance energies (E = 3.75MeV corresponding to E Uh = 105 MeV for 28 Si + 28 Si). 
At a short distance of the folding potential, i.e., with highly overlapping densities, 
DDM3Y gives weakly attractive potential. At the normal density, for example, the 
density-dependent factor f(E,p) reduces the interaction strength by a factor about 
3/4, compared with the original g(E,r), while it is enhanced by a factor 1.2 at the 
half density, i.e., at the contact region. 

The folding-model potential, however, is considered to be accurate only in the 
tail region of the nucleus-nucleus interaction. In the region where nuclear-density 
overlap goes beyond the normal density, it is considered to be not accurate enough. 
Hence, in addition to the folding potential with the nucleon-nucleon interaction, 
we introduce a phenomenological repulsive potential, which would originate from 
the effects of the Pauli principle among nucleons belonging to the interacting nuclei 
respectively, or from compression effects due to the overlapping density. We estimate 
strength of the repulsive potential due to the compression of nuclear density, from 
the equation of state of nuclear matter, i.e., from the binding energy as a function of 
nuclear density. One may think that the picture of the density overlap is doubtful in 
low energy, but the folding model does not take into account density redistribution, 
so it is consistent to account higher densities in the overlapping region. Anyhow, 
what we are interested in is the dynamics of two interacting nuclei in high spins where 
strong centrifugal forces dominate. Therefore, the long-range part of interactions is 
crucially important, but not the short-range part, which is treated more or less in a 
phenomenological way. 
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The folding potential is defined as usual, 

U(R) = J dr 1 J dr 2 pi(r 1 )p2(r 2 )v(r 1 2), 

r 12 = R + r 2 -rt, (2-24) 

where R is the relative vector between the interacting nuclei and T{ are referred to the 
centers of the nuclei, respectively. The long-range attractive part of the interaction 
potential in the molecular frame I4ttr is obtained from U(R) by taking the vector R 
to be parallel to the z'-axis and by taking orientations of the density distributions 
of the constituent nuclei with respect to the molecular frame. By using Fourier 
transformation, 

^ttr = ^2^r'y /m (fi) J dkk 2 jt(kR) J dkl^ m (k)v(k)^{k)p2(-k), (2-25) 
Im 

v(k) = J dre ikr v(r) = 4vr J drr 2 j (kr)v(r), (2-26) 
Pi(fe) = I dr'e ikr ' Pi (r'). (2-27) 



The density distribution pi(r') in the molecular frame is related to that in the body- 
fixed frame, i.e., to that in the principal axes of the constituent nucleus; by Euler 
rotations, Pi{r') = pf{r'[) = TZ(aij3i^i) pf '{r') , where pf{r'() is the density distribu- 
tion in the principal axes and therefore 

pfir'l) = Yl fnW)Y n ffl) (2-28) 

l=even 

with the assumed axial symmetry of each constituent nucleus. So the Fourier trans- 
form joj(fe) is given with the Euler angles included as parameters, 

Pi(k) =Y / i l Mk)J2 D m'o(^iPai)Yi m '(k), (2-29) 

/ m' 

Pl(k) =4ir J dr^jtikr)^). (2-30) 

Inserting Eq. ()2-29j) with i = 1 and 2 into Eq. ()2-25j) . we obtain the final form of 
the interaction potential as a function of the internal variables (R, a, /3x, (5 2 ) in the 
following, 



Kttr(i?,a,/3i,/3 2 ) = W 3 i l '- l "- l l'i"(l'l"00 I 10) 

l'l"l 

xF vvn (R)G vvn {a,l3i,l3 2 ), 
F n »i{R) = j dkk 2 ji(kR)v(k)p v {k)pi»(k^ (2-31) 

G vm (a,p u 02) = Et" 1 )"^ 2 " <W(/Tm - m | 10) 

m>0 

x cos(2ma)c4 (/3i)cC (/? 2 ). 
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It should be mentioned here that ji does not appear in the final expression due 
to the D-function with one magnetic quantum number being zero which originates 
from the axially-symmetric density distribution in Eq. (|2-28|) . and that a\ and «2 
are combined into 2a = a± — Q2 due to the fact that the vector R is parallel to 
z'-axis, i.e., the magnetic quantum number associated with R is zero. The Coulomb 
interaction is also folded, together with nuclear interaction v(E,p,r) of Eq. (|2-2ip . 

We assume the density profile of pf{r'[) to be the Fermi distribution with 
pf(r'l) = p /[l + exp{(r" - i? N (r'/))/a N }], #N( r i) denoting the radius of the de- 
formed nucleus. As for the deformation of the constituent 28 Si nuclei, the existence 
of the hexadecapole deformation (04 = 0.18±0.02) is suggested from coupled-channel 
analyses for the elastic and inelastic neutron scatteringPB Therefore, we take the 
radius of each nucleus as R^ir") = tqA^ [1 + /3q12o(^') + P^Xioi^i)] including two 
parameters /3q and /3h for the deformations, the values of which are determined to 
be —0.46 and 0.22, respectively, according to the suggested value for the ratio /3q//3h 
and their magnitudes adjusted with the B(E2) value of the ground-rotational band of 
28 Si.E3 The value of r is taken to be 1.03fm from the textbook of Bohr-Mottelsonj23 
and on to be 0.48fm to reproduce the RMS radius of the ground state. 

Next, we proceed to the effect of density overlap in the inner region, where the 
folding potential is not expected to be adequate. An overlapping of the densities 
brings about a higher nuclear density than the normal one, which gives rise to a 
binding energy loss of the interacting system in addition to the attractive folding 
potential. We take into account the effect as a repulsive potential to be added to 
the folding one given in Eq. (|2-24|) . The volume with higher density depends on 
the configurations of the constituent nuclei, especially on their relative distance. 
Actually, the overlapping of two nuclei produces nuclear density from zero to twice 
of the normal density. An accurate calculation of the effect, therefore, is rather 
laborious. We propose a simple approximate way. If we assume the density profile 
to be of sharp cut-off or with a very small diffuseness, an overlapping volume has 
always twice of the normal density. So the short-range repulsive effect is expected to 
be proportional to the overlapping volume, and it would be simulated by a potential 



where the primes on the densities indicate Fermi distributions with a small diffuseness 
ap. The strength Vp of V rep is chosen in the following, referring to the Equation of 
State (EOS) of nuclear matter. Thus the total interaction potential is given by 



where V^ttr denotes the usual folding potential defined in Eq. ()2-24p . The repulsive 
potential looks like a folding potential of the zero-range interaction, but has the 
primed densities instead of the normal density distributions. Of course, we can 
utilize a merit of the form of Eq. (|2-32p in the actual calculations. 

To determine the strength Vp, we use EOS of the nuclear matter, i.e., a binding 
energy loss per nucleon As for twice of the normal density which is calculated under 




(2-32) 



Vint — ^attr + *z 



(2-33) 
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j=o 




10 11 



R(fm) 



the condition of complete overlap at the 
R = limit. Without Coulomb en- 
ergy the value of Ae can be taken to 
be 7 ~ HMeV^S from the values of 
the nuclear compression modulus = 
180 ~ 240MeVp3 which is suggested by 
the experiments on giant monopole res- 
onances. Hence the values ap = 0.25fm 
and Vp = 330MeVfm 3 are obtained to 
reproduce Ae = 9MeV in the 28 Si + 28 Si 
system. Radial forms of the folding po- 
tential are shown in Fig. 3, for the sta- 
ble geometrical configurations (parallel 
equator-equator ones, see the next sec- 
tion), where the effective potentials for 
J = and J = 38 are displayed. Details 
of the folding potentials about their de- 
pendences on ap and Vp are already dis- 
cussed in Ref. 120]) . where the effects of hexadecapole deformation in 28 Si nuclei are 
also investigated. 

§3. Dinuclear structures of the 28 Si+ 28 Si system 

3.1. Stable configuration of the oblate-oblate system with high spins 

In order to know dynamical aspects of multi-dimensional internal motion, we 
calculate the effective potential with specified spin J and K, defined as follows: 



Fig. 3. The radial forms of the effective 
potentials in the parallel equator-equator 
configuration of 28 Si + 28 Si, V JK (R) = 
Vi n t(R, tt/2, tt/2, 7r/2)+T r ' ot (J, K) for spins 
J = and J — 38 with K = are shown. 



v JK = Mnt (/J, a, h) + r r ' ot ( J, K) . 



(3-1) 



R(fm) 




12 L 1 1 1 1 ~l 

10 20 30 40 50 60 70 80 90 100110120 
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Fig. 4. The effective potential energy Vjk for the 8 Si + Si system with J = 38 and K = 
is displayed, for the R — = /?a) degrees at a — tt/2. A local energy minimum exists at 
R — 7.6fm. Contours are in MeV. The figure is essentially the same as the energy contour map 
in Refs.[!l} andl2TJ|). 
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In Fig. 4, an R — fafa = fa) energy surface, i.e., Vjk(R, n/2, /3, /3) is displayed for 
J = 38 and IT = 0. We find a local minimum point at fa = fa = n/2 and R = 7.6fm, 
namely, at the equator-equator(E-E) configuration, with a rather deep potential well 
around the equilibrium. We mention that the some expressions, numerical results 
and figures in this section are already published in Refs. [19]) and I20p. but we show 
them for explanation. 

In Fig. 5(a), the a-dependence of 
Vjk in the E-E configuration at the 
equilibrium distance is shown. (Note 
that our definition for the domain of the 
variables is < a < n and < fa , fa < 
n.) We find that the a-dependence is ex- 
tremely weak. Another point is that we 
have two local minima at a = and n/2. 
Those two configurations are, however, 
exactly the same, namely, parallel E-E 
configuration (z"-axes of the constituent 
nuclei are parallel). Therefore it is nec- 
essary to impose symmetry on the wave 
functions. In Fig. 5(b), /3-dependences 
of Vjk with fa = fa are compared be- 
tween at a = and at a = 7r/2, where 
solid line is for a = n/2 (the cross sec- 
tion of Fig. 4 at R = R c = 7.6fm) and 
dashed line for a = 0. (Note that con- 
figurations with fa = fa 7^ 7r/2 at a = 




50 



70 90 110 
(3 (degree) 



130 



Fig 



5. (a) the a-dependence of the effective 
potential Vjk with J = 38 and K = 0, for 
the E-E configuration at R — R e = 7.6fm. 
(b) Vjk versus /3 = /3i = @2 at a — and 
at a — it/2, which are displayed by dashed 
and solid lines, respectively. The figure is 
the same as published in Refs. 119(1 and!20p. 



are not the same as those with the same 

/3j-values at a = n/2, but are the same as those with fa = tt — fa at a = tt/2.) The 
fa well at a = n/2 is seen to be rather shallow, compared with that at a = 0. Hence, 
despite the weak a-dependence of Vjk in the E-E configuration, we have significantly 
a-dependent restoring force for /3-motions around the E-E configuration. 



3.2. Harmonic approximation and normal modes with a specified K 

In order to solve normal modes for four variables (R, a, fa, fa), we expand Vjk 
into a quadratic form for R, fa and fa, at the equilibrium E-E configuration, while 
for a we keep its dependence exactly in terms of cos(2ma) series, such as those given 
in the interaction potential of Eq. (|2-3ip . Then the effective potential is expressed 
as 

V JK (R,a,fa,fa) = Vj K (Re,a,~) + ^(R- R*) 2 

+ hl l {a)Af3l + hf{a)Af3l 
+ kf{a)AfaAfa + {higher order), (3-2) 
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where Afc denotes — tt/2. k % l{a) denotes the second derivative d 2 VjK/dj3idj3j, 

k l p{a) being equal to k 22 (a). Although k l ^{a) is a coefficient of A(3iA(3j in the 
expansion, it is a function of a, i.e., we take into account a-dependence of the 
coefficient, in addition to the a-dependence of Vjx(Re, ct, |, j). As k\ (a) consists 
of cos(2ma) series with m = even including zero, the major part is a constant ko from 
m = 0. We write fci (a) = fc« 2 (a) = ko + ^2(0), ^2(0) being a sum of contributions 
from terms with m = even > 0. 

We introduce new coordinates in order to eliminate cross products of /?i and 
02 both in Ti n t and in the quadratic expansion of Vjk- The new variables describe 
butterfly and anti-butterfly modes as follows: 

/3+ =(Afa + A(3 2 )/V2 = (ft + /3 2 - tt)/\/2, 

/3_= (Ap 1 -Ap 2 )/V2 = (p 1 -p 2 )/V2. (3-3) 

Furthermore the inertia masses of three variables a, f3 + and /?_ are approximated by 
the values given at the E-E configuration. Combining the kinetic energy operator 
and the expanded effective potential, the total hamiltonian is given as follows: 

H = H Q +T' C + {higher order), (3-4) 



H = H R + H lmgl (l3 + ,/3-,a), (3-5) 

^-tm + T*-^ (3-6) 
# ang i(/3 + ,/3_,a) = #+(/?+, a) +#_(/?_, a) 

h 2 d 2 TA 7T 7T. /n 

~4J^2 + ^ ( e,Q! '2'2 ) ' (3 " 7) 

FT ffi ^- ^Z 1 1 1±cos2q \ ^ fc±(q) 2 , . 



where + or — sign of ± in Eq. (|3-8p corresponds to the /?+ and /3_ degrees of freedom, 
respectively, with k+(a) = ko + £2(0:) + k 12 (a) and fc_(a) = &o + ^2 (a) — k^ 2 (a). 

Now we solve the Schrodinger equation with the hamiltonian Hq for the internal 
four degrees of freedom, which is separated into two parts. One is the hamiltonian 
Hr for the radial motion and nothing but that of a simple one dimensional harmonic 
oscillator. Another is -ffangi for the angle variables a, /3+ and /3_, which is also almost 
separable into H + of H- of /?_ and the remaining hamiltonian for a. H + and i?_ 
again represent harmonic oscillators, although the masses and the restoring forces 
depend on a. Hence we analytically obtain wave functions for H± and their energy 
quanta hw± with the frequencies 



,1 l ± cos 2a \ 

,jJ - r - ~ \ihJn)[j+ —2 ) . (3-9) 
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Taking into account those vibrational energies from the /3-degrees of freedom, we 
introduce a reduced potential for the a-motion, and obtain the Schrodinger equation 
for the a-motion as follows: 



- I - 5 ^ + V JK {R^a,-,-) + EP n+>n _ 



a 



4>(a) = E angl (, 



a 



(3-10) 



where En + ,n-{ot) denotes vibrational energy + l/2)?ki; + + (n_ + l/2)hu^ from 
H + + H^, added as a part of the reduced potential. Note that, in order to obtain 
analytic form of hia± in cos(2ma) series, we expand square root in Eq. (|3-9p supposing 
loq = ^Jk${l/I + 1/ fiR 2 ) to be the leading term. Accordingly, we consider a solution 
4>(a) of Eq. (|3-l(jp to be described by cosine and sine functions of a, i.e., Fourier 
series, as the reduced potential Vjx(Re,C(, § , § ) + En +:n _(a) is described by a sum 
of cos(2ma). Then Eq. (|3TU|) is reduced to a secular equation, which is easily solved. 
Thus the eigenenergy of the system is given as follows, specified by the quantum 
numbers (n,n + ,n_, K, (u, ir a )), 



n 2 



E J (n,n + ,n_,K,(u,7r a )) =E (R C ) + — 



J( J + 1) - K 2 - 1 K 2 - 2 



fiR 2 
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+ (n + + n- + l)Hu + E2(TT a ), (3-11) 

where v denotes a dominant frequency of the a-motion with 7r Q for the parity con- 
cerning the reflection at the equilibrium of a = ir/2. The first and second terms in 
the r.h.s. of Eq. ()3Tip are constant energies from the interaction potential and the 
centrifugal energy included in Vjk at the equilibrium, respectively. (n + + n_ + 1)^0 
and£°(7r a ) are the vibrational energies for the /3-motions without the a-dependence 
and the energy for the a-motion, respectively. 

There is a selection rule K ± v = 
even for the a-motion. Because of the 
parity and boson symmetries, n+ can 
be taken to be larger than or equal to 
n_. For the /3-vibrational modes, we 
have a rule (-1)"++"- = (-1)^ due 
to the symmetry of each constituent nu- 
cleus under the space inversion. Details 
of the symmetries of the molecular sys- 
tem, the wave functions and the selec- 
tion rule are given in Appendix B. The 
resultant states are summarized in Ta- 
ble I. Note that the eigenfunction of the a-motion is not necessarily the internal 
rotation specified with a single i^-value, and that mixing over allowed ^-states is 
expected. 

In Fig. 6, molecular normal modes of 28 Si + 28 Si with spin 38 are displayed, 
classified with the /C-quantum numbers. The twisting-mode excitations associated 



Table I. Molecular states allowed by the selec- 
tion rule, specified by the /("-quantum num- 
ber, the /3-vibrational quanta (n+, n_) and 
v for the a-motion. 



K 


(n+,n_) 


V 





(0,0), (2,2) 


0, 4, 8,... 





(2,0), (4,0), (4,2) 


0,2,4,6,... 


2 


(0,0), (1,1), (2,2) 


2, 6,10,... 


4 


(0,0), (1,1), (2,2) 


0, 4, 8,... 


2, 4 


(2,0), (4,0), (4,2) 


0,2,4,6, ... 


1, 3 


(b0), (2,1) 


1, 3, 5,... 
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55 
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28 Si + 28 Si J = 38 



(2,1) 



(2,1) 



(1.1) 



(2,0) 



(2,0) 



(2,0) 



(1,0) 



£=4 (1,0) _ 
(0,0) 

®(B)®(B)®(B)®(B)©(B) 

^K=0^ ^K="H ^K=2 * *- K=3^ fc K=4 * 




butterfly 




anti-butterfly 



Fig. 6. Molecular normal modes for the Si + Si system for J = 38. The quantum states are 
specified by (n,n+,n-,K,(i>,ir a )), where n = is given except for one level(n = l,v ~ 0) 
displayed with dashed line. The quanta (n+,n_) of the /3-motions are given below the levels, 
and K at the bottom, (t) and (b) marks assigned in the lower part of the figure indicate the 
twisting rotational mode and the butterfly modes, respectively. Also given above some levels 
with K — are dominant values of the quantum number v for the a-motion. On the right-hand 
side, butterfly and anti-butterfly motions are illustrated. The figure on the l.h.s. is the same as 
published in Refs.Ql} and 120)1 . 



with the a-degree are obtained, and indicated by marked (t) at the bottom. Also 
given above each level with K = and (t) is the dominant quantum number of v for 
the a-motion, which means the a-motion is approximately described by a single term 
cosz^a. The butterfly and anti-butterfly vibrational modes are indicated by marked 
(b). A pair of quanta (n + ,n_) is given below the levels. All those are due to the 
internal degrees of freedom, i.e., intrinsic excitations. Apparently the i^-excitation 
and the twisting rotational mode appear to be lower than the /3-vibrational modes. 
The excitation energy for K = 2 is very small, smaller than lMeV, and even those 
for K = 4 or v = 4 are smaller than 3MeV. 

In Fig. 7(a), a few examples of wave functions for the a-motion are exhibited, 
where the /3-modes are in the zero-point oscillation (dashed line) or the 2-quanta 
excitation of butterfly (solid line). We see that, with zero quanta for the /3-modes, 
the amplitude is wriggling around the value of the unit, the equilibria a = and 7r/2 
being slightly favored. (With exact v = we have a constant behavior. Weak v = 4 
mixing exists.) With 2 quanta for the butterfly mode, however, we find surprisingly 
strong concentration around the equilibrium of a = tt/2. In Fig. 7(b), we inspect 
the reduced a-potential for quanta (2,0). Compared with the potential for (0,0), 
we find that the minimum at a = disappears, and the potential well at a = tt/2 
is extended to wider region, which sustains the localization of the amplitude. One 
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may wonder why the difference between 
a = and tt/2 exists. The reason is as 
follows: at a = 0, due to the definition 
of f3±, /3-motion with (n + ,ra_) = (2,0) 
does not imply butterfly excitation but 
anti-butterfly one with 2 quanta. Such 
a characteristic of the f3± coordinates 
gives larger excitation energy for (2, 0) 
at a = than at a = tt/2. In Fig. 7(c), 
the energy quanta ftw± versus a are 
shown, where we are able to confirm the 
point. Returning back to the dinuclear 
configuration, for a configuration with 
Pi = fa < tt/2, for example, we ob- 
tain a butterfly one at a = tt/2, such 
as displayed in Fig. 2, while at a = 
we obtain an anti-butterfly one with the 
same values of /3j. Hence the localization 
around a = tt/2, seen in Fig. 7(a), indi- 
cates nothing but a realization of a phys- 
ical butterfly excitation. Thus, we are 
able to classify the levels in Fig. 6 into 
two groups, i.e., the twisting mode and 
the butterfly (or anti-butterfly) mode, 
respectively. Some examples of the wave 
functions for the normal modes are ex- 
plicitly given in Appendix C. 




3 0.5 1 1.5 2 2.5 3 



a 

Fig. 7. (a) Wave functions for the a- 
motion for J — 38 and K = 0. Those 
with the zero-point oscillation (0, 0) and 
the butterfly excitation (2, 0) for the /3- 
degrees of freedom are displayed, re- 
spectively, (b) The reduced a-potential 
Vjk(R c , a, n/2, tt/2) + E^ +>n _ (a), (c) a- 
dependences of the /3-energy quanta hu±. 
The figure is the same as published in 
Refs.ni]) and [HI- 



§4. Rotational motion at extremely high spins with triaxial deformation 



One of the characteristic features of 
the spectrum obtained theoretically is a 
series of low-energy i^-rotational exci- 
tation due to axial asymmetry around 
molecular z-axis, which is in contrast 
with the 24 Mg + 24 Mg caseP' 1 ^ One 
can understand the reason immedi- 
ately from Fig. 8, where the upper 
configuration( 24 Mg + 24 Mg) has axial 
symmetry as a total system, but the 
lower one for 28 Si + 28 Si has axial asym- 
metry. Thus K is not a good quantum 
number, namely, we expect the eigen- 
states are K-mixed. 

A triaxial system preferentially ro- 




Y 



Fig. 8. Equilibrium configurations of two di- 
nuclear systems. The upper portion is 
for 24 Mg — 24 Mg and the lower one for 
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tates around the axis with the largest moment of inertia. By the definition of the axes 
in the lower panel of Fig. 8, we have the moments of inertia as Ix > Iy » Iz, due to 
the nuclear shape. Thus the system, which is seen as two pancake-like objects( 28 Si's) 
touching side-by-side, rotates around X-axis normal to the reaction plane. Such a 
motion is called as wobbling. 

We extend our molecular model so as to include couplings between states with 
different i^-quantum numbers. As a result, we will obtain new low- lying states due 
to the triaxial shape of the equilibrium configuration. The Coriolis terms in the 
molecular hamiltonian bring those couplings. However, in practice, we do not treat 
the Coriolis terms explicitly, but we diagonalize the hamiltonian of the asymmetric 
rotator to obtain the rotational spectrum. 

The Coriolis terms in the molecular hamiltonian in Eq. (|2-16p gives an impression 
that those are quite different from the asymmetric rotator. So the effect of the 
Coriolis coupling terms will be examined later in §4.2, to show that the molecular 
hamiltonian reduces to the asymmetric rotator hamiltonian in the sticking limit. 

4.1. Analyses by asymmetric rotator 

We describe the rotational motions of two pancake-like objects( 28 Si's) touching 
side- by-side by means of the asymmetric rotator. Generally its hamiltonian is written 
as follows, with the moments of inertia about the intrinsic axes I x , I y , and I z , 
respectively; 



trot 




(44) 

where J x , J y and J z denote the components of the angular momentum operator 
along the intrinsic axes of the body-fixed frame (the same operators as defined in 
Eq. (|2Tip ). 7 av , A and Ik in Eq. (|4-2p are related to I x , I y and I z by 

f-|(f4). <«> 

± a,v a \ A x ± y / 
11111 r* ^ 

A = -T + - = T--> (44) 

1 x J av 1 y 1 av 



1 1 1 

IK Iz lav 



(4-5) 



By using lowering and raising operators, J+ and J_ of the angular momentum in 
the body-fixed frame, we obtain 



[ rot 



where J± = J x ±iJ y , respectively, which give rise to couplings between different i^'s. 
The coupling strength is given by the coefficient 1/A, which is proportional to the 
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difference between 1/I y and 1/I X - In an intuitive understanding, the rotation around 
x-axis is lower in energy than the rotation around y-axis due to I x > I y . When the 
energy difference between the rotations around the molecular x- and y-axes is larger 
than ET-excitation energies, the iC-mixing is expected to be rather large. In other 
words, an energetically favored motion, i.e., rotation around x-axis would be realized 
by the if-mixing. 

In order to obtain an accurate description of this triaxial rotator, as it is well 
known for polyatomic molecules, we diagonalize the hamiltonian with an inertia ten- 
sor of the axial asymmetry, which gives rise to mixings of ET-projections of the total 
spin jJ3S) The resultant motion should be called as "wobbling mode" PD The energy 
spectrum is displayed in Fig. 9(b), compared with the spectrum without K- mixing in 
Fig. 9(a). Now the states of low lying iT-series are not the eigenstates by themselves, 
but are recomposed into new states. It is very interesting that we again obtain sev- 
eral states including the K = component as a result of iT-mixing, which should 
show up themselves in the scattering. Those states are closely located in energy and 
so in good agreement with several fine peaks observed in the experiment. It should 
be noted here that due to the lack of the enough information about the deformations 
of the total system, we assumed the same parameter for the coupling strength A for 
the molecular ground-band states and for the butterfly states, although the extent of 
asymmetry is generally different in each band. As for the magnitudes of the moments 
of inertia (I x , I y , I z ), we estimated them as follows. We assumed a constant value 
for the relative distance, i.e., for jiR 2 , and adopted R = R e = 7.6fm. For the con- 
tributions from the moments of inertia of the constituent 28 Si nuclei, we estimated 
them about the y- and z-axes from the excitation energy of their 2f state. As for 
the contribution about x-axis, we assumed a factor 4/3 larger than those about the 
other axes, due to the distribution of the nuclear density of 28 Si. For calculations 
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Fig. 9. Energy spectra of the 28 Si + 28 Si system for J = 38. (a) Molecular normal modes without 
AT-mixing. (b) After A'-mixing, with indications of the modes under the levels. 
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about the moments of inertia of dinuclear systems, details are given in Appendix D. 

As an analytical prescription, in the high spin limit (K/J ~ 0), the diagonaliza- 
tion in the i^-space is found to be equivalent to solving a differential equation of the 
harmonic oscillator with parameters given by the moments of inertia. Thereby, the 
solution is a gaussian, or a gaussian multiplied by an Hermite polynomial, 
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where the second term on the r.h.s. is due to the coupling energy between the states 
with AK = 2, which is approximated by K/J = 0. The third term is due to the 
energy of the harmonic oscillator, with the energy quantum, 
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Its excitation should be with n = even due to the symmetry between the \K > and 
| — K > components. Note that by the approximation J(J + 1) ~ J 2 for the second 
term of the r.h.s. of Eq. (|4-9f) . the first and second terms can be amalgamated into 
J(J + l)h 2 /2I X , which reminds that the moment of inertia of the rotation is I x . 



Considering I x ~ I v , i.e., I 
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the energy quantum fwj of Eq. ()4-10l) 



is equivalent to that of the wobbling formula given in Ref. 35). 

In order to calculate angular correlations we use those analytic forms in Eq. (|4-7|) , 
which is simple and intuitive way to understand the extent of .ff-mixing. Of course 
we can utilize numerical values obtained in the diagonalization procedure, but the 
values are almost the same as those given by the analytic form. For the lowest state 
fo(K) of Eq. (|4-7p . we have the wave function for the wobbling ground state as 
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where in general, xk can be any molecular mode of triaxial deformations, such as 
the ground-state configuration (parallel equator-equator one), the butterfly mode 
and the anti-butterfly mode. The magnitude of b estimated by Eq. (|4-8p is 1.85, for 
example, for the values of the moments of inertia used in the calculations for the 
energy spectrum in Fig. 9. This is the largest value expected, because we assumed 
a static configuration there, in which the zero-point motions of the twisting and 
butterfly modes are neglected. A note for the wobbling wave functions is given in 
the last part of Appendix C. 
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In Fig. 10, theoretical energy levels of the Si + Si system are compared with 
the experimental data in the resonance energy region.^® From the left, (a) shows 
the molecular ground band followed with the wobbling excited states, (b) an excited 
band due to the twisting motion with K = 0, (c) the butterfly mode with wobbling 
and (d) the anti-butterfly mode with wobbling. Molecular configurations are well 
stable by the barrier up to J = 40, while with J = 42, an existence of the molecular 
resonance state is unlikely, as the zero-point energy of the radial motion is over the 
barrier top. Levels with J = 38 are connected by thin lines for eye-guide. On the 
right-hand side, the experimental data are displayed, where J = 36 ~ 40 indicate 
the spin assignments for the broad bumps P We see the density of the resonance 
states in the data is well reproduced by the calculated eigenstates, which are due to 
the wobbling motion and the excitations of the internal modes, such as butterfly etc. 
Note that the existence of the excited states of the wobbling motion as resonance 
states depends upon the stability of the triaxial structure, and hence the numbers 
of those excited states taken up in Fig. 10 are not definitive. 

From J = 36 up to J = 40, the anti-butterfly mode appears higher than the but- 
terfly one as is discussed in §3.2. For those anti-butterfly states we do not display the 
excited states of the wobbling motion, because each eigenstate of the anti-butterfly 
mode appears as the excited state of the butterfly state in the a-motion, and thus the 
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Fig. 10. Dinuclear spectrum of the 28 Si + 28 Si system in the resonance energy region is displayed. 
Levels with J = 38 are connected by thin lines for eye-guide. From the left, the resonance levels 
theoretically obtained, where (a) the members of the molecular ground band and of the wobbling 
excited bands, (b) an excited band due to the twisting motion with K — 0, (c) the butterfly 
mode with wobbling and (d) the anti-butterfly mode with wobbling. On the right-hand side, 
the experimental data are displayed, where J = 36 — 40 indicate the spin assignments for the 
broad bumps.® We selected resonance levels from the narrow peaks in the elastic and inelastic 
excitation functions, according to a statistical analysis on their correlation, with an indicative 
bar for the available energy region of the dataP** 
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configuration of the anti-butterfly mode is not enough triaxial. On the other hand, 
with J = 34 the anti-butterfly mode is lower in energy than the butterfly one, and 
thus we display the excited states of the wobbling motion for this mode. The reason 
of the lower excitation of the anti-butterfly mode is as follows. With relatively-low 
angular momentum, the constrain by the energy well around the equator-equator 
configuration becomes rather weak, and the stability of this configuration is not well 
guaranteed. We found that the equator-equator configuration is not at the local 
energy minimum below with J = 32. For example, with J = 30, the stable con- 
figuration for the molecular ground state is an antibutterfly-like one of (3 ~ 60°, 
the equilibrium distance R c of which is much smaller than that with J = 34. Such 
softening of the energy surface occurs with J = 34, which gives rise to lowering of 
the anti-butterfly mode. 

4.2. Comparison between the molecular-model hamiltonian and the asymmetric 
rotator's 

The wobbling motion associated with "X-mixing" is an important aspect of 
the rotation of the asymmetrically deformed nucleus in high spins. In order to 
investigate such an aspect, we have introduced the asymmetric rotator in addition to 
the molecular model, because the rotator model is simple to understand the essential 
feature of rotational motions. As the rotator model is based on more or less rigid 
intrinsic structure, it is interesting to know how the simple rotator model is related 
to the molecular model. 

A triaxial system preferentially rotates around the axis with the largest moment 
of inertia. By the definition of the axes in the lower panel of Fig. 8, we have the 
moments of inertia of the total system as Ix > W » Iz due to the configuration. 
Thus the total system, which is seen as two pancake-like objects touching side- 
by-side, rotates around the X-axis which is normal to the reaction plane. In this 
context, the magnitudes of the moments of inertia is crucially important; the large 
contributions to Ix from the third moments of inertia I3 of the two constituent nuclei 
are expected. So here we study two examples of the molecular model, one of which 
is with I3 = 0, and the other is with ^3 / 0. 

We take up a resonant system consisting of " a spherical nucleus and a deformed 
nucleus" . In order to see only the rotational motion, we assume that the two nuclei 
are bound and stay at a constant relative distance R, which reduces the degrees 
of freedom of the system. For the first example (case 1), an axially-symmetric 
deformation is assumed for the constituent nucleus, in which .Z3 is taken to be zero. 
Thus the coordinates are taken as (g^) = (9i, $2, #3, /3), where 9i denote the Euler 
angles for the rotation of the molecular axes. The internal degree of freedom is 
described with (3. Those conditions are taken to be corresponding with the degrees 
of freedom in §2. On the other hand, in the second example (case 2), an axially- 
asymmetric deformation is assumed, i.e., I3 is not zero, which gives rise to a degree 
of freedom 7. And then the coordinates are (%) = (#1, 62, #3, 7). Of course, the 
introduction of the 7-degree of freedom is an extension from the description in §2, 
which is expected from the nuclear density distribution illustrated in Fig. 8. 

Typical configurations are illustrated in Fig. 11, where in (a), (b) and (c) the 
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Fig. 11. Some examples for the geometrical configurations of the system consisting of a spherical 
nucleus and an axially-symmetric deformed nucleus. The molecular z'-axis is set to be parallel 
to the z-axis, as the whole rotation by f2(9i, #2, #3) will start with those configurations. In 
(a), the symmetry axis of the deformed nucleus (z"-axis) is parallel to the z'-axis, due to 
the Euler rotations with (a,/3) = (0,0). In (b), the z"-axis is parallel to the x'-axis, due to 
(a, /3) — (0, 7r/2). In (c), the z"-axis is parallel to the y'-axis, due to (a,f3) — (tt/2,w/2). 

configurations are set with (a, /3) = (0,0), (0, 7r/2) and (7r/2,7r/2), respectively. 
The configuration (a) is axially symmetric for case 1, while those of (b) and (c) 
are axially asymmetric (the same shape). Note that the shape of the constituent 
deformed nucleus is prolate, but the figure is useful both for the prolate nucleus and 
for the oblate nucleus, of course. We describe the orientations of the principal axes 
of the constituent deformed nucleus with the Euler angles (a, /3), in which the degree 
of freedom associated with a is essentially the same as that associated with #3 of 
the total system, and thus a does not appear in (g$). Actually, the rotations of the 
whole system and the constituent nucleus are described by the relation, 

Q n (aJ^) = <(a,/?, 7 )r2 M (0i, 2 ,0 3 ), (4-12) 

where Q n and Q' n denote Euler rotations for the constituent deformed nucleus with 
respective angles, and Qm denotes rotations of the molecular axes. On the r.h.s. 
of Eq. (|4-12p . Q' n denote the successive rotation after Qm\ firstly the axes of the 
constituent deformed nucleus rotate up to the directions of the molecular axes by 
Qm-, and secondly they rotate referring to the molecular axes by Q' n . Since the 
successive rotation is decomposed into Q' n {^)Q' n {j3)Q' n {a), we have the rotations 
with angles #3 and a around the same axes obtained after ($l> $2)' Thus we can 
take the angle of the third rotation of the whole system simply to be 63 + a, which 
involve the freedom a. Here we regard the value of a as the initial condition for the 
starting configurations before rotation, which are displayed in Fig. 11. 

We write the classical kinetic energy with the angular velocities, and then we 
replace these angular velocities with time derivatives of those coordinates (%). The 
classical kinetic energy is expressed in the form T = | Yl 9ijQiQj an d we quantize it 
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by using the general formula for the curve-linear coordinate system. The quantum 
mechanical expression for the kinetic energy is given by 

where g and g -1 denote the determinant and the inverse matrix of (gij), respectively. 
The metric tensor (gij) is composed with the submatrices g ro t for the whole rotational 
degrees and gi n t for the internal degrees of freedom as 

<*> - ( 'To Z ) • (4 - 14> 

where gc denotes the nondiagonal part which corresponds to the Coriolis coupling, 
t gc being the transpose of gc- The way to obtain those components of g^ is described 
in detail in Ref. I18p . Here we briefly see their definitions: 

Srot = V(0 2 , 03)(I„ + P, i)InR(a, /?, j))V(9 2 , 9 3 ), (4-15) 
9mt = *n^, 7 )J n y(/3, 7 ), (4-16) 
<?c = V(0 a , 3 )*i2(ar, /3, 7 )I»^GS, t), (4-17) 

where V(#2) $3) an d V(^,7) denote the transformation matrices between the deriva- 
tives of the Euler angles and the angular velocities of the molecular axes and of 
the constituent deformed nucleus, respectively, R(a, f3, 7) being the rotation matrix. 
Ifj, denotes the inertia tensor for the two constituent nuclei as point-masses, i.e., 
the diagonal moments for x- and y-axes being fiR 2 , while I n denotes the inertia 
tensor for the constituent deformed nucleus in its principal axes, respectively. The 
moments of inertia of the constituent deformed nucleus (the diagonal elements of 
I n ) are taken as follows: for the axially symmetric nucleus (case 1), I\ = I2 = Io 
and ^3 = 0, and for the axially asymmetric nucleus (case 2), I\ = I2 = Iq and the 
value of ^3 being not zero. For the latter case, in general, I\ 7^ I2 may be used 
for the static asymmetric deformation, as was tried in the theory of the asymmetric 
rotatorpSJbut for simplicity we avoid this tedious calculations. Our consideration is 
focused on the appearance of the 7-degree of freedom associated with ^3 / 0, and 
for this purpose the assumption I\ = I2 brings no problem. Note that 7 is spurious 
for the deformed nucleus with the axial symmetry, i.e., for case 1 with ^3 = 0. With 
the aid of mathematical software, we can easily obtain the elements g and (g~ l )ij, 
for example, g = (^R 2 ) 2 Iq sin 2 j3 sin 2 O2 for case 1, and g = (//i? 2 ) 2 /o^3 sin 2 (3 sin 2 62 
for case 2, respectively. 

As the classical kinetic energy consists of three parts, i.e., the rotation of the 
whole system, the internal motions and their couplings, the quantum mechanical 
operator for the kinetic energy T is also given as a sum of three terms, 

f = f TOt + f int + f c . (4-18) 

Naturally the term T rot is associated with the rotational variables (0±, 62,63), T{ nt 
with the internal variables (/?, 7) and Tq with both. According to the derivation, T ro t 
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is expressed by the partial differential operators of 6%. We combine those differential 
operators into angular momentum operators J[. Thus we obtain the following kinetic 
energy terms, 

f ** = Y E Hjj'A ( 4 - 19 ) 

l<i<3 
l<i<3 

where the matrix [i is the submatrix given later, and J[ are the angular momentum 
operators in terms of the Euler angles of the molecular frame, as usual, which is 
already given in Eq. ()2-lip . The coefficients fjj>i j 3X6 determined due to the moments 
of inertia corresponding to the geometrical configuration, and are given in terms of 
the parameter a and the internal variable f3 as follows; 



(4-20) 



1 

Note that the above expressions of are the same between two examples, which 
are obtained under the assumption 1\ = I2 = I®. 

The internal kinetic energy operator is associated with the variable f3 for case 1. 
We obtain 

= -2 U + w ) w - t U + w ) + v ' (4 ' 21) 

where the second term on the r.h.s. is the additional potential due to the new volume 
element dV = d/3 instead of the original dV = fiR 2 1 sin /3df3. For the asymmetrically 
deformed constituent nucleus with 1% ^ (case 2), we have an additional term 
associated with the 7-degree of freedom, i.e., 

h 2 ( 1 / 1 1 \ 1 1 r) 2 

lUVn) = JUOT - y{-^(^™ 2 ^ W ) + T 3 }w (4 ' 22) 

where indications (/?) and (/3, 7) on Tj nt are for distinction between cases 1 and 2. 

The Coriolis coupling operator Tq consists of coupling operators between the 
variables (^1,^2,^3) and j3 for case 1, i.e., 

TciOn P) = j - sina( - i— ) J[ + cosa( - 4 J, (4-23) 

where the derivative operators of 9i are rewritten with the angular momentum oper- 
ators J[. As for case 2, the couplings are between the variables (#1, #2, 63) and (/3, 7), 
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for which we obtain 

fcVf.M) = fb(»,;« - <*»«( - i£) i; + 8inQ ( _ ;|.);<} 

-^'"(rtX-^H (4 ' 24) 

Now we investigate those kinetic energy operators obtained. Firstly we proceed 
with case 1. Generally the hamiltonian of the system is composed of the kinetic 
energy and the interaction between the constituent nuclei V(/3), i.e., 

H = f mt + {f^ + V(P)} + f c . (4-25) 

Without the coupling Tq, the internal motion is determined by the eigenvalue equa- 
tion, 

H int xifi) = {T^t + V{p)} X iP) = E X {P). (4-26) 

For the case of a very strong confinement by the interaction V(/3), the motion as- 
sociated with the /3-degree is expected to approximately follow Eq. (|4-26p . And we 
write the rotational motion of the system, including To, as 

H rot = f mt + f c . (4-27) 

To analyze the rotational motion given by H rot , it is important to assume a domi- 
nant configuration such as the equator-equator one of the 28 Si + 28 Si system, which 
appears in low energy due to the interaction between the constituent nuclei. For we 
have all the possible rotational motions generally to appear in the energy spectrum 
of .Hrot j which make the problem extremely complicated. The dominant configu- 
ration assumed in this investigation is an axially asymmetric one with the strong 
confinement which is seen in Fig. 11(b) and/or Fig. 11(c). Thus we put ft ~ ir/2 for 
fiij of T ro t in Eq. (|4-20p . i.e., cot/3 to be zero. As for a, its value is not essential, 
and we can choose any value for a. The Coriolis coupling Tc(^;/3) of Eq. (|4-23p 
reduces conveniently into one term, for the configuration with a = illustrated in 
Fig. 11(b), and we obtain 

H 2 f J( J + 1) ( 1 1 \ / d 2 \ 2 / ti . 

where the first term of the r.h.s. is the angular momentum J 2 replaced by J( J + 1), 
and the last term gives the Coriolis coupling. Note that for another configuration 
with a = tt/2 illustrated in Fig. 11(c), — J[ appears instead of J' 2 . In general, the 
Coriolis coupling induces X-mixing associated with the vibrational excitations, which 
is exemplified by the third term of the r.h.s. of Eq. (|4-28P ; by using the creation and 
annihilation operators, a* and a of the vibrational motion obtained by Eq. (|4-26l) . 
the term is rewritten as ~ (a* — a)(J_ — J+), where = J[ ± J' 2 are the lowering 
and raising operators of the iT-quantum numbers. This description corresponds to 
that of the molecular hamiltonian given in §2. 
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On the other hand, there is a picture in which the internal degree f3 is frozen due 
to strong adhesion between the constituent nuclei, i.e., the sticking limitW^ In order 
to obtain relation to the rotator hamiltonian, we approach from such a picture, in 
which the sharing of the angular momenta is determined classically, i.e., 

J = (^R 2 + I )cj, L = liR 2 lj, S = I u, (4-29) 

where u denotes the angular velocity of the whole system. Corresponding to the 
above relations, quantum mechanical ones are derived in Eq. (|AT8p of Appendix A 
as follows, 

S/Iq = L/fiR 2 = J/ifiR 2 + J ), (4-30) 

where S, L and J denote the spin of the constituent deformed nucleus, the orbital 
and the total angular momenta, which are defined as the rotations on the same 
plane. Although the relations in Eq. (|4-30p are limited among the rotations around 
the same axis, they are approximately applicable to the present analysis, because 
we are considering the rotation of the whole system around the largest moment of 
inertia; for example, the axis of the rotation is y' for the configuration in Fig. 11(b). 
Therefore S is associated with the /3-degree, and J corresponds to J' 2 with respect 
to the configuration in Fig. 11(b). Following Eq. ()4-30p . we use a relation, 

.d_ 

l d(3J 

and rewrite the total kinetic energy T = T rot + T; n t + Tc • Then the result turns out 
to be the hamiltonian of the asymmetric rotator, i.e., 



-) /h = J' 2 /{iiR 2 + I Q ), (4-31) 



2 \ iiR? fiR 2 + I I 

Note that due to the configuration in Fig. 11(b), we put a = and f3 = ir/2 for 
Hij of r rot in Eq. (|4-2U|) . and Tc(^i] P) of Eq. (|4-23|) . respectively. Note also that we 
drop the additional potential in T- int {j3) (the second term of the r.h.s. of Eq. (|4-2ip ). 
because the new volume element is not applied for the spin degree of freedom (3. 
Thus we successfully obtained the asymmetric rotator hamiltonian. However, T in 
Eq. (|4-32p is too simple, because the /3-degree is frozen and disappears. So this way 
does not exactly correspond to what is shown in §4.1. There we have the hamiltonian 
for the activated /3-degrees in addition to the rotator one. Although the above T 
in Eq. (|4-32p is well correspondent to the kinetic energy operator given in §2.2 and 
gives a rotator hamiltonian, eventually, the results of the diagonalization of Tq in 
§2.2 may not give rise to the energy spectrum given in Fig. 9. 

Next we examine the case 2, in which two internal degrees of freedom (/?, 7) are 
treated. In Fig. 12, the configuration with the oblate deformed constituent nucleus 
is illustrated corresponding to Fig. 11(b), where a larger value of ^3 than that of 
Iq is assumed. So the axis with the largest moment of inertia is x' in this case. 
According to the configuration, we again put a = and f3 = ir/2 into T rot , Ti nt (/3, 7) 
and Tc(#i; /3, 7), respectively, which gives 

A in N ^ f ( 1 1 \ ( 92 l \ I 1 1 \ 5 2 1 , lnns 

IUCA7) = -y{ ( To + ^) ( w + ~ 2 ) + (5 + W ) 7? \, (4-33) 
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and 

^.;A7) = -|,{(-4)4 + (-^)i;}. (4-34, 

We again proceed in the picture of the sticking limit. Relations similar to Eq. (|4-29j) 
are assumed as follows: 

J = {^R 2 + I 3 )u, L = fj,R 2 u, S = I 3 u, (4-35) 

where u denotes the angular velocity of the whole system again, the rotation being 
around the x'-axis. Due to the same relation as Eq. (|4-30|) except Iq being replaced 
with I3, we use a relation, 

{-%) lh = J'xl^R 2 + h). (4-36) 

By Eq. (|4-36p . only the 7-degree is frozen out, and as a result we obtain the total 
kinetic energy operator, 

f & r (jjf , (4) 2 , & \ 

2 \ixR? + Iz ixR? I j 
h 2 ( l l \ / d 2 i\ h 2 ( . d \ f/ 

Thus we finally obtain the asymmetric 
rotator hamiltonian accompanied by the 
vibrational mode (3. It is noted that 
the equilibrium of the /^-vibration is as- 
sumed to be (3 = tt/2 in the derivation. 
In this model, the pole orientation of 
the constituent deformed nucleus is fluc- 
tuating around the direction of x'-axis 
seen in Fig. 12, while the whole system 
with the moment of inertia fiR 2 + 13 ro- 
tates approximately around the x'-axis 
with rather confined configuration. This 
is just the same picture adopted in §4.1, 
where the two constituent nuclei keep 
in touch and are sticking in the rota- 
tional motions along the same axis as 
the whole rotating system. As for the 
moments of inertia of the system, strong 
confinements due to the nucleus-nucleus 
interaction are supposed to give rise to 
induced deformations and /or a neck for- 
mation of the constituent nuclei, which bring nonzero moments J3 of the constituent 
nuclei in addition to the original moments Jo of the deformed nuclei with the axial 
symmetry. Thus, we expect that the above analysis is physically meaningful. 



Z, Z' 
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Fig. 12. An example for the geometrical con- 
figurations of the system consisting of a 
spherical nucleus and a nucleus with oblate 
deformation, for intuitive understanding. 
The molecular z'-axis is set to be paral- 
lel to the z-axis, and the symmetry axis of 
the deformed nucleus (V'-axis) is parallel 
to the x'-axis, due to (a, 0) = (0, 7r/2). 



Resonances in Si + Si. I 



31 



§5. Summary 

The interaction between two nuclei is described with the internal collective vari- 
ables, i.e., the orientations of the poles of the constituent nuclei in the rotating 
molecular frame. In the dinuclear system with oblate-deformed constituent nuclei, 
an equator-equator touching configuration with the parallel principal axes is found 
to be the equilibrium at high spins. In the 28 Si + 28 Si system, the relative distance 
between the two 28 Si nuclei is 7 — 8fm, indicating a nuclear compound system with 
hyperdeformation. The barrier position is 9 — lOfm, greatly outside from that of 
usual optical potentials. Molecular configurations are well stable by the barrier up 
to J = 40, while with J = 42, an existence of the molecular resonance state with 
narrow widths is unlikely, as the zero-point energy of the radial motion is over the 
barrier top. This theoretical maximum spin is in accord with the bumps observed 
in grazing angular momenta. 

Couplings among various molecular configurations are taken into account by 
the method of normal mode around the equilibrium configuration, which gives rise 
to the molecular modes of excitation, such as the radial vibration, the butterfly 
motion, the anti-butterfly motion and so on. The twisting mode {y = 4) is found 
to be the lowest excitation. Vibrational energy quanta for the butterfly and the 
anti-butterfly modes are about 4MeV, but the excitation energies of those modes 
have to be twice, 8MeV, since states of K = even with one vibrational quantum are 
not allowed due to the boson symmetry. Thus, the energies are close to those for the 
radial excitation. Although the excited state of the radial mode is not bound in the 
present calculations, the possibility of the radial-mode resonance is not completely 
excluded, because it is likely that the interaction between two 28 Si would be more 
attractive than the present folding potential with the frozen density. 

A triaxial system preferentially rotates around the axis with the largest moment 
of inertia. By the definition of the axes in the lower panel of Fig. 8, we have the 
moments of inertia of the total system as Ix > W » Iz- Thus the total system, 
which is seen as two pancake-like objects touching side-by-side, rotates around the 
X-axis which is normal to the reaction plane. As the axial symmetry is slightly 
broken, wobbling motion appears in that way. 

We extend our molecular model so as to include couplings between states with 
different X-quantum numbers. Usually, the Coriolis coupling terms are diagonalized, 
but we do not treat them explicitly. In practice, we use the asymmetric rotator as an 
intuitive model. By the diagonalization of the rotator hamiltonian in the X-space, we 
obtain new low-lying states due to a triaxial shape of the equilibrium configuration. 
In the high spin limit {K/J ~ 0), the diagonalization is found to be equivalent 
to solving a differential equation of the harmonic oscillator with spring constants 
given by the moments of inertia. Thereby, the analytic solution is obtained to be a 
gaussian, or a gaussian multiplied by an Hermite polynomial, which is a useful tool 
for the analyses of the molecular states with the triaxial configuration. 

Since the Coriolis terms in the molecular hamiltonian appear to be quite different 
from the asymmetric rotator, it is necessary and meaningful to study the relations 
between the molecular hamiltonian and the asymmetric rotator's. The analysis turns 
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out that the hamiltonian of the molecular model with the 7-degree of freedom reduces 
to that of an asymmetric rotator in the sticking limit. 23 ' Thus the intuitive use of 
the asymmetric rotator is warranted, and it provides a very simple understanding 
with easy calculations of the effects of if-mixing on the energy spectrum. 

Finally it should be mentioned that an extension of the molecular model is 
possible so as to include the 7j-degrees of freedom. For example, possible 7- vibrations 
of the constituent nuclei could be taken into account. But we do not include those 
surface vibrations and the corresponding 7j-degrees of freedom, considering that the 
dominances of the members of the ground-state band in the decays are reported for 
the 28 Si + 28 Si system® and for the 24 Mg+ 24 Mg system,^ respectively. A molecular 
model with two asymmetric rotators of the constituent nuclei is not pursued for 
the moment, which does not appear rewarding for elaboration. Furthermore, the 
constituent nuclei are expected to be strongly confined to form the whole deformed 
system, in which the 7j-degrees of freedom are approximately frozen. Hence, we 
adopt the asymmetric rotator for the whole system as a sticking limit of the 74- 
degrees of freedom. 

We have intuitively expected that the moment of inertia Ix is the largest for 
the configuration in Fig. 8. Namely, we implicitly assume I\ 7^ I2 and ^3 7^ for the 
moments of inertia of the constituent nuclei in their principal axes, due to additional 
deformations likely induced by the interactions at the contact configurations of the 
two nuclei, while in §2, without the induced deformations, we have assumed the 
axial symmetry of 28 Si and the intrinsic moment I3 = 0. (The latter gives ly > 
Ix >> Iz) The dynamical process of the transition between those two states of the 
constituent nuclei with ^3 = and ^3 7^ is an interesting problem, which should be 
clarified in future. The nuclear structure with large induced deformations may be 
close to that obtained by the Nilsson-Strutinsky calculations for 56 NiP2J However 
the experiments exhibit the nature of the dinuclear complex in resonances, i.e., the 
dominance of binary decays,® which suggests the contact of the two 28 Si nuclei is 
not violent enough for rearrangements of the nuclear structure of the constituent 
nuclei in the molecular model. 

The molecular states obtained in the present paper are expected to be the origin 
of a large number of resonances observed, and hence theoretical analyses have been 
made for the angular distributions and the angular correlations. The results have 
been compared with the recent experiment performed in Strasbourg®'^ to give 
good agreements with the data, which will be given in the succeeding paperP^ 
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Appendix A 

Relation between the Coordinate Systems and the Angular Momentum 

Operators in the Molecular Model 

In the present appendix, the de- 
scription of the angular momenta for 
the total kinetic energy operator in the 
molecular model is ascertained with re- 
spect to the coordinate system. For this 
purpose, we take up a two-body problem 
in which one body is deformed and has 
own rotational degrees of freedom. In 
the laboratory frame the description of 
the kinetic energy of the system is rather 
simple as we consider the energies from 
the rotational motion of the two-body 
relative vector (orbital motion) and the 
spin degrees separately. However in the 
body- fixed frame, i.e., in the molecular 
frame, the description is not so easy, be- 
cause the total system is not simple as 
a rigid rotator. Firstly, the coordinates 
for the molecular frame and those of the 
internal degrees of freedom associated with the frame should be chosen appropri- 
ately. Secondly, although we can obtain the kinetic energy operator T by using the 
formula for the curve- linear coordinates which is given later in Eq. (|A 3|) . we have 
three parts of T, i.e., the rotational energy, the energy associated with the internal 
degrees of freedom and the couplings between them, the roles of which have to be 
clarified. As is shown below, the property of the operator, for example, to be the 
total angular momentum or to be the orbital angular momentum, is determined 
not only by the coordinates for the rotational motions themselves but also by the 
moments of inertia associated with. This is natural as the angular momentum in 
the classical description. Two examples are given; one is for the present molecular 
frame, and the other is for a new molecular frame, in which the couplings between 
the whole rotation and the internal motions disappear. 

Consider a resonant system consisting of "a spherical nucleus and an axially- 
symmetric deformed nucleus" . In order to see the rotational motion, we assume that 
the two nuclei are bound and stay at a constant relative distance R, which reduces 
the degrees of freedom of the system. To limit the degrees of freedom to be "two 
dimensional" , we again assume that the symmetry axis of the deformed nucleus is 
in the reaction plane. As is illustrated in Fig. 13, referring to the space-fixed axes, 
the orientation of the relative vector of the two constituent nuclei is described by the 
angle 0$, an d the angle of the orientation of the symmetry axis is denoted as oq. The 
classical kinetic energy is given with two angular velocities ui = 6q and u) n = do, as 
T = 1/2 • (^R 2 lo 2 + In^n), where \x and I n denote the reduced mass of the two nuclei 



y 

A X" 




Fig. 13. The Coordinates of the system 
which consists of a spherical nucleus and 
a deformed nucleus with axial symmetry. 
Both the relative vector and the symmetry 
axis of the deformed nucleus are assumed 
to be parallel to the xy-plaxie. 
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and the moment of inertia of the deformed nucleus, respectively. The corresponding 
kinetic energy operator is given by 

h 2 h 2 

■ " l 2 + ^tS\ (A-i) 



2[iR 2 21, 

where L and S denote the orbital angular momentum and the spin of the deformed 
nucleus, as usual. Note that by definition, the angular momenta are those for the 
one dimensional rotations, i.e., they are given by L = —id/dOo and S = —id/dao, 
respectively. 

In the description by the molecular model, we take the molecular x'-axis which 
is parallel to the relative vector of the two nuclei, as illustrated in Fig. 13, and the 
coordinate is denoted as 9 (the angle 9 is the same as 9q, i.e., 9 = 9q). And the 
angle of the orientation of the symmetry axis is described referring to the molecular 
x'-axis as a = Qo — 9q. Then due to oj' n = u n — oj with cu' n = a and lu = 9, the 
classical kinetic energy is written as 



T= l - 



\ (/total + 2I n uu/ n + I n oj' n 2 ), (A-2) 

where /total denotes the moment of inertia of the whole system given by /total = 
liR 2 + I n . Replacing these angular velocities with time derivatives of the coordinates 
(<?«) = we write a classical kinetic energy in the form ^ ^2gijqi<ij- And then 

we quantize it by using the general formula for the curve-linear coordinate system, 

t = -^y J- A ./sor 1 ) (a-3) 

1 j 

where g and g^ 1 denote the determinant and the inverse matrix of (gij), respectively. 
In this case, the metric tensor is given by 

M = ( " R ) + In (A-4) 



and hence the inverse matrix is obtained as 

(a-Y l -( l/ ^ R2 _1/ ^ 2 ^ fA-5) 
{g v } -{ -1/^R 2 l/^ + l/I n J' (A5) 

As the classical kinetic energy consists of the three parts, i.e., the rotation of the 
whole system, the internal motions and their couplings, the quantum mechanical 
operator for the kinetic energy T is also given as a sum of three terms, which appears 
to be 



f = M J- A _ J d l_ + (J_ + (A . 6) 

2 \ fiR 2 d9 2 fiR 2 89da ^ V ^R 2 ^ ij da 2 J 1 ; 

We define the total angular momentum Jg and the operator for the internal (ro- 
tational or vibrational) motion of the deformed nucleus S a , by Jg = —id/ 89 and 
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S a = —id/da, respectively, and rewrite Eq. (|A-6j) with those operators as follows, 

ti 2 

T = —{jlliiR 2 - 2J e S a /fiR 2 + (1/^R 2 + l/I n )S 2 a } (A-7) 



h 2 „ , 2 



2 



(j e -S Q ) 2 + — S 2 ,, (A- 



which corresponds to Eq. (|A-lj) and exhibits the role of the Coriolis terms in the 
molecular model. Note that the coordinate 9 for the molecular x'-axis is the same 
angle as 9q, but the role of 9 is different from 9q. The molecular axis x'{9) represents 
the motion of the whole system with the moment of inertia /total > while 9q is the angle 
of the relative vector between the two constituent nuclei, which represents the orbital 
motion. The process by using the formula Eq. (|A-3[) of the general quantization for 
the curve-linear coordinates gives a simple example for the kinetic energy operator 
described in §2, which clarifies the property of the molecular coordinate 9. 

We can introduce the angular momentum operators Jq and S a , of course, by the 
direct transformation for the differential operators. We obtain, due to the relations 
between the arguments (0q = 9,ao = 9 + a), 

d 89q d dao d d d 

89 = ~Wd9~o~ + ~d9~dc7 a = W + da^' ( 

d_ = d9o_d_ + dao_d_ = _d_ (A-10) 

da da 89q da dao dao ' 

and accordingly we can confirm the relations, Jq = L+S and S a = S. As for the wave 
functions, let us start those with the eigenvalues M and m for L and S, respectively. 
The total wave function \P is defined by x l r (9Q,ao) = Ne lM s ° e %ma ° , the arguments 
of which could be replaced by the coordinate transformation as (9q = 9, «o = & + 
a), and accordingly we obtain &(9,a) = Ne iMe e im{ - e+a ^ = Ne i ( M+m ^e ima . Note 
that the function e vma is the same one obtained from the operation of the unitary 
transformation of the whole rotation R n (9) = e ~ ie ( n 'S) G n the spin function in 
the laboratory system, i.e., R n (9)e ima ° = e lTna , where n denotes the unit vector 
normal to the plane. (For the general rotations in three dimensional space, the 
transformations are described by D-functions.) The resultant part e l ( M + m ) s properly 
corresponds to the wave function for the degree of freedom of "the whole rotation 
of the system" , and thus we again confirm that the eigenvalue J of Jg satisfies the 
usual rule J = M + m. 

In the molecular model, firstly we consider the whole rotating system with a 
stable (equilibrium) configuration expected, and secondly we investigate the internal 
degrees of freedom associated with it. With the strong nucleus-nucleus interaction, 
the motions of the constituent nuclei may be perfectly confined, and hence we some- 
times consider the internal degrees of freedom to be frozen, i.e., the sticking UmitW^ 
Thus it is worth while looking the kinetic energy of the molecular model in the 
sticking limit. In the classical kinetic energy Eq. (|A-2p . we put oj' n = (uj = u n ) 
and accordingly we obtain the energy of the rotator T = I tota io; 2 /2. However, in 
the quantum mechanical expression in Eq. (|A-6|) and/or (|A-8p . it is clear above that 
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the molecular model does not directly correspond to the sticking model, because the 
rotational energy is not given by J 2 fr? /2I tota \. To make a model corresponding to the 
sticking model in quantum mechanics, we define a new coordinate 0, "the angle for 
the center of the moments of inertia" in stead of the Euler angle for the molecular 
z'-axis, as = (fiR 2 ■ 6q + I n ■ ao)/(/j,R 2 + J„). The other coordinate is again "the 
internal angle" a' = ao — 9q, which is the same as in the molecular model, and the 
corresponding moment of inertia is given by /internal = I^R 2 ■ In/i^R 2 + In)- The set 
of the coordinates (0, a') gives the classical kinetic energy expression without the 
Coriolis coupling term, as 

T = » (-^total^ 2 + -^internal^ 2 ), (ATI) 

where Q denotes the angular velocity 0. Thus we obtain the kinetic energy operator 



2-^total ~ 2/jnternal 

where by putting S' a = we reach the kinetic energy of the rigid-rotator type. 
Note that the choice of those coordinates follows the usage of the center of mass 
coordinate and the relative vector for two-body problem. Unfortunately this set of 
the coordinates would be limited on the rotations in a plain, because it is not easy 
to define "the center of the moments of inertia" for the multi-dimensional internal 
rotations. Now, by the direct transformation for the differential operators, we again 
calculate operators J© and S' due to the relations 



which appear as 



#o = - In/(vR 2 + In) ■ a', (A- 13) 

ao = + fiR 2 /(fiR 2 + I n )-a', (A- 14) 



Odd , . _ N 



30 89q dao 

d I n d iiR 2 d 



da 1 fiR 2 + I n dd ' fiR 2 + I n da ' 

Thus again we have a usual description for the total angular momentum associated 
with as a sum of the orbital angular momentum and the spin, i.e., Jq = L + S. 
As for the spin for the internal rotation S' a , the definition turns out to be 

S' a = -id/da' = (fiR 2 -S-I n - L)/{ixR 2 + I n ). (A-17) 

For the sticking limit, we put S' a = 0, and then we obtain 

S/I n = L/fxR 2 = J/ifiR 2 + I n ), (A- 18) 

which is known as the rule of angular momentum sharing in the sticking modelPS 
Thus the coordinate system taken up here properly gives the whole rotation and the 
internal motion without the coupling as a quantum mechanical description of the 
sticking model. 
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Appendix B 

Symmetries of the System and Construction of the Wave Functions 

Here we deal with symmetry properties of the system and their associated re- 
strictions on the wave function, from which selection rules for quantum numbers are 
deduced. Following those results, practical expressions for wave functions are given 
in the next Appendix C. 

Firstly we note the coordinate transformation rules for boson and parity opera- 
tions in the molecular frame. Here we do not describe how to obtain the rules. One 
could refer the derivations given in the Appendix D of Ref. I18|) . 

Boson symmetry 

We have the exchange operator Vi 2 , which acts on both the molecular coordi- 
nates and the internal variables, and transforms them as follows; 



V l2 : {6i,62,e 3: a,R,(3 u (3 2 )-*(TT + 6i,tt - e 2 ,-0 3 ,a,R,7r - p 2 ,n ~ Pi). (B-l) 



V : (0i, 2 , 03, a, R, /3i,/3 2 )— K?r + 0i, vr - 2 , vr - 3 , -a, R, (B-2) 



Wave functions of the system with good symmetries 

Since the axial symmetry of constituent nuclei is assumed, the variables 7» are not 
necessary. Each nucleus has positive parity, and thus its density profile is invariant 
under space inversion. Accordingly, the basis wave function D J MK (6i)xK{R, ct, Pi,P 2 ) 
should be invariant under the inversion operation upon a constituent nucleus, 



Before we examine the symmetry for Zj, it should be noted that the transformations 
Zj affect the rotational variable 03 as well as the internal variables a, Pi and p 2 , 
because orientation of the molecular x'-axis changes according to a change of the 
orientation of a constituent nucleus. For example, we take up Z 2 : (ai,/?i) — >(&i + 
27r,/?i), which should be equal to unity, because it gives just 2tt rotation of one 
constituent nucleus around the molecular z'-axis. By operating X 2 on 03 and a of 
D J MK (9i)xK{R,oi, Pi, /3 2 ), according to 3 = (ai + a 2 )/2 and a = (ai - a 2 )/2, we 
obtain a cyclic condition including the factor from the transformation on 03, 



D J MK (9 i ) XK (R,a,Pi,p 2 ) = {-l) K D J MK {e i ) XK {R,a + Ty,Pup2)- (B-4) 



Now in order to examine symmetries about the inversion operations Zj, we set 
trial wave functions concretely. By introducing harmonic approximation with the 
variables p+ = (A0x + Ap 2 )/^/2 and /3_ = (APi - Ap 2 )/V2 with Ap { = ft - vr/2, 
the internal motions are described with 



Inversion symmetry (parity) 

The inversion operator V acts as follows; 




(B-3) 
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where n, n + and n_ denote oscillator quantum numbers, respectively. For 
XK(R,a,fa,fa) of Eq. (|B-5j) . Eq. (|B-4[) means as 

0jc(a) = (-l)*>jf(a + ir). (B-6) 

Note that a-dependences in y?+ and originate from those in the oscillator hamil- 
tonians H±{f3±, a) for the /3+ and /3_ degrees of freedom, which are almost separable 
from the a-degree but not completely, as is seen in Eqs. (|3-5|) — f|3-8j) in §3. Natu- 
rally, the periodical property of the a-degree of freedom is of period ir. Furthermore 
due to the geometrical identification of the configurations specified with a = n/2 and 
a = 0, we have a relation between butterfly function (£>+ + (/3 + ,a) and anti-butterfly 
one (p~_(J3-,a). According to Eq. (|3-9p . we have those oscillator energies Huj + = 

h^k + {a){l/I + (1 + cos2a)/fiR 2 } and Huj^ = Hy/k-(a){l/I + (1 - cos 2a) /^R 2 }, 
respectively, where k + (a) and k-(a) denote spring moduli for respective modes. The 
moduli are defined by the coefficients of AfiiAfij in the harmonic expansion, and have 
been written in the text as k + {a) = ko + k 2 {a) + k^ 2 (a) and k-(a) = ko + k 2 (a) — 
k l2 {a), respectively, where ko is a constant, and /02(a) and k l2 (a) consist of cos(2ma) 
series with m = even and m = odd, respectively. Since cos 2 (a + vr/2) = — cos 2a 
and this is also the case in kjj 2 (a), i.e., k^ 2 (a + vr/2) = —kV 2 (a), we have relations 
fthj + (a + 7r/2) = faj_(a) and ftw-(a + 7r/2) = faj+(a). Thus the oscillator functions 
,(/?+, a) and <p~_(/3-,a) also follow the same relations, such as 

<p+(P,a) = <p- l (p,a + Tr/2) = <p+(p,a + <ir), (B-7) 

where the variable /3 denotes /3 + or /3_. As for the transformation on A fa, the oper- 
ator Xi is the inversion by the definition, and so we have X\ : (#3, a, j3 + ,fa) — >-(#3 + 
7r/2,a + vr/2,-/3_,-/3 + ) and Z 2 : (0 3 , a, P-)-±{6 z + vr/2, a - vr/2, /?_ , /3+). By 
utilizing the above relations, we perform symmetrization of the basis wave function 
about Xi. Starting with the wave function D{ 4K (6 i )xK(R, a, fa, fa) with Eq. flBJ)} 
for xk(R,o:, fa, fa), we obtain 

«P A = (1 +X 1 ){1 +X 2 ) • I)' Jf (e i )/„(l«)^(a) V +.03 +J a)^-_C9_,a) 

= £>^(^)/n(^){l + (-l) n++n - X }^n + n„(a,/3+,/3-), (B-8) 

with 

h K n+n-(a,P + ,fa) = (f) K (a>)<Pn + (P+,aL)iPn-(P-i a ) 

+(-i) K ct )K {a + tt/2)^+_ (/?+, a)<^ + (/?_, a). (B-9) 

From the phase in the braces of Eq. (|B-8p . we obtain the selection rule for /3-mode 
quanta as 

(-1)"++"- = (-1)* (B-10) 

Note that the second term of the r.h.s. of Eq. (|B-9p originates from X\ and X 2 , with 
n+ and n_ exchanged. 
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Subsequently, we perform boson symmetrization and parity-projection by oper- 
ating ^(1+7^12) and l) pr P} to We obtain the symmetrized wave functions 
as follows: 



^ ~ D J MK (0i)f n (R)\h Kn+n _ (a, p + ,p_) + (-iy +K h Kn+n _ (-a, -?+,?-) 
+{ _ 1)P+ j-k d .J m _ k{9i)Ur) 

x \l Kn+n A-^(3 + ,(3-) + (-lf +K h*_ Kn+n _(a,-(3 + ,l3.)\. (B-ll) 

Note that as the V12- (or V-) operation gives —K for the total rotation, we set —K in 
Eq. ()B-lip for /i^ n+n _ (a, /?+, /3_)-functions consistently, as follows. Since Eq. (|B-6[) 
indicates the sign for the cycle for \K\, we define as (a) = ^jj-(a). With the real 
functions ^xipt), the effects — K in /i_R- n+n „ (a, f3 + , /?_ )-functions appear only in the 
phase (— i) K of the second terms of the r.h.s. of Eq. (jB-9[) for K = odd as complex 
conjugate. We further reduce hxn+n- (a, /?+, /?-) in Eq. (|B-llj) with the arguments 
—a and/or — /?+. Since the reduced potential in Eq. (|3-10|) has reflection symmetries 
at a = and a = tt/2, we are able to classify 4>k(o) by parities with respect to those 
points, i.e., 

4>K{-Ct) = 7T a=0 • 4>k{oc), 

cp K (tt/2 - a) = vr a=7r/2 • X (tt/2 + a) , (B- 12) 

where 7r a= o and vr Q , =7r / 2 denote the parities with respect to the reflection points, 
respectively. Note that a relation 7r a=7r /2 = (— l) x vr Q , = o is known because of 4>k(ci + 
7r) = (—1) K 4>K(a) by Eq. ()B 6|) . Due to a) = fvjj{a) for each f3 + - or /3_-mode, 
we also k now <p+,(J},—a) = <£+,(/?, a) and y~,(/3, -a) = <p~,(fi,a). By applying 
Eq. (|B-12j) to hxn+n- ipt, /?+,/?_) of Eq. (|B-9|) . we rewrite the internal wave functions 
in Eq. (|B-lip . For the first line, for example, we obtain 

h K n + n- («, P+,P-) + (-l) P+K h K n + n. ("«, -/?+,/?-) 

= {1 + n a=0 ■ (-lf + "-}{^(a)^+ + (/3 + ,a)^_ (/?_,«) 

+ {-i) K <t> K {oi + 7r/2)<p+_ (J3+, a)<p- + , a)} , (B-13) 

where the rule Eq. (jB-10jl . (— 1) = (— \^ n ++ n - \ s used. Hence, relations are obtained 
as ir a= o ■ (— l)P+ n - = l ; and for even parity states we have 

7r a= Q = (-l) n ~, 

7r Q=7r / 2 = (-l) n+ , (B-14) 

which specify parities of a-mode in connection with the /3-mode quanta. Under the 
parity selection rule in the a- motion in Eq. (|B-14j) . we can rewrite the functions of 
Eq. (|B-13p into hicn + n- («> P+, P-)- Thus, the final form of the total wave function 
with the symmetries is as follows: 

*x ~ D J MK (9 i )f n (R)h Kn+n ^{a,(3 + ,^) 

+{-l) J + n +D{, K (9 i )f n {R)h. Kn+n _ (a, p+, £_), (B-15) 
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with the definition of hjc n+n _ (a, /3+, in Eq. ()B-9p and with the selection rules 
given in Eq. (|B40l) and Eq. 

Restrictions on quantum numbers 

As a summary, we have selection rules given in Eqs. (|B-10|) and (|B-14p . some 
conditions for the wave functions such as given in Eq. (jB-7[) and the total wave 
function given in Eq. ([Bd^]) with Eq. (fB~9]) . In the following, we give some practical 
selection rules reduced from those relations. 

Due to the cyclic condition <f>Kip- + tt) = (— ^-) K 4>k{ol), if we expand 4>K(ct) 
with periodic functions such as cosi^a or sinra, we have a general restriction for 
rotational quantum numbers (K ± v) = 2m, m being an integer. And then, with 
a specified K, one of the two parities of Eq. (|B-14p is enough for specifying 4>k((x) 
to be cosine type or sine one, the other selection rule being automatically fulfilled. 
Because of the symmetry of each constituent nucleus under space inversion, n+ 
can be taken to be larger than or equal to n_. In the case that n + is equal to 
n_, we have K = even from (-1)™++"- = {-1) K of Eq. (|B40jh And Eq. (TEH))) 
turns out to be {^(a) + (— i) <Pk{oi + ir /2)}tp^, (f3 + , a)ip~,((3- , a) , which suggests 
{K ± v) = 4m. (See also Eq. (|C-2j) for if-, i/-rules.) As for K = 0, we obtain 
the phase rule (— l)^ +n + = 1 from Eq. (jB- 15|) . i.e., n+ = even and ra_ = even for 
J = even. With if = v = 0, the state has no a-dependence and hence the second 
term of Eq. (|B-lip become to be the same as the first term, which gives the phase 
rule (— l) p+J = 1, i.e., J = even for the positive-parity states and J = odd for the 
negative-parity states, respectively. 

Appendix C 

Explicit Expressions of Wave functions for the Normal Modes 

In the present Appendix, we give some examples of explicit expressions for the 
wave functions. As is shown in Appendix B, the internal wave functions are approxi- 
mately a product of f n {R), (a), ip^ (f3 + ,a) and p~_(/3-,a), which are essentially 
oscillator wave functions except for 4>k{ol)- The property of a-motion may be ro- 
tational, or may be vibrational, depending on the strength of the reduced potential 
which confines the alpha-degree of freedom. Note that the reduced potential is deter- 
mined by the interaction between two nuclear surfaces as well as by quantum states 
of the other degrees of freedom. Especially in the molecular ground state, where 
the additional potentials from the normal-mode motions such as the butterfly mode 
is weak, the property of 4>k(o<-) is determined by the interaction. The confinement 
potential obtained from the folding model is weak as is shown in §3, but its reality is 
not confirmed yet. Hence we present both types for 4>k{ol) functions, i.e., cosine and 
sine series as well as gaussian functions. First we take up the rotational type and 
adopt cosine series for <^(a), assuming K, v = even and n + ,n_ = even. Note that 
e lva may be convenient for describing /i^ n+n _ (q, /3 + , /3_) in Eq. ()B-9p . but it does 
not fulfill Eqs. (pd2l) and (jB~T4)l . Note also that for K = odd, h Kn+n _ (a, /3+, 
includes both cosine and sine functions due to the shift by tt/2 in 0^(a). 

Now we write down the total wave function according to Eqs. (|B-9p and (jB-15p . 
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For positive-parity states with K = even and n + = even (n_ = even, 7r Q=7r /2 = 
vr a= o = 1 ), by applying cos v[a + 7r/2) = (—l) u / 2 cos z^a for 4>k{®- + tt/2), we have 

^ ~ [^(^) + (-l) J ^ j _ K (0 i )]/ n (i?)^n + n_(«,/3+^-), (C-l) 

with 

h Kn+ n-(a,(3+,(3-) = ^2 C ^ cosva { < i Pt + {P+^ a ) L PnAP-^ a ) 

u=even 

+(-l)^+J/3+,a)^ + (/?_,a)}. (C-2) 

If we take a single value of u, the a-motion is, of course, rotational. The quantum 
state (n,n + ,n_, K,u) = (0,0,0,0,0) corresponds to the molecular ground state. 
Actually the solutions obtained in §3 is close to that with v = 0. For K = 0, the 
next rotational state is with v = 4, namely, the twisting rotational mode. The explicit 
functions for /i^ n+n _ (a, /3_) of those states are given as 

v = : /ioooK/3+,/3-) = y^c^o (/?+, a) (/?o (/?_,«), 

z^ = 4: /i 00 o(a,/3+,/3_) = y^cos4a^(/3 + ,a)^o (/3_,a). (C-3) 

Due to the weak a-dependence of the folding potential, solutions receive small mix- 
ings over v. The molecular ground state has the v = 4 component about 3% in 
the probability, as well as the v = component in the v = 4 excited state. Those 
coefficients CVs are adopted for the calculations of the partial decay widths and 
the angular correlations, with the approximation of the constant vibrational quanta 
for hp + and Hp_ = ^/3_) 5 the value of which is taken to be 4MeV. As for the 

relative motion, it is completely separated from the (a, /3+, /3_)-degrees of freedom, 
and it is described by oscillator wave function f n {R) with the center at the equilib- 
rium distance R e . According to the experimental resonance energy E cm = 55.8MeV 
for J = 38, the molecular ground state with n = is a suitable assignment as the 
theoretical eigenenergy is 51.5MeV. With n = 1, the radially-excited state appears 
at higher than experimental energy, and it may not be observed in experiments, be- 
cause of the broad decay widths expected. On the other hand, there is a possibility 
that a stronger interaction appears due to induced deformations, which give rise to 
lowering the eigenenergy of the radially-excited state to 55.8MeV. Note that char- 
acteristics of probability distributions among the decay channels receive no effect 
from the choice of the radial motion with given n, which is completely separated 
from the angle degrees of freedom. Only the magnitudes of the partial decay widths 
commonly become larger as we take a higher n-value. 

If we coherently sum up over v in Eq. (|C-2p . the a-motion can have a localized 
property such as a zero-point oscillation. For example, for simplicity, we take up 
n + = n_ = 0, and then the selection rule for K and v is K ± v = 4m with m being 
an integer, because two terms in the braces of the r.h.s. of Eq. (|C-2p must have the 
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same sign. We can write simply as 

h K oo(a,P+,P-) = ^2 C„ cos vatp^(/3 + ,a)ip Q3_,a). 

K+v=im 

Now, if we adopt gaussian type coefficients for C v by 



C v 



1 \l/4 
1 + 8v0 ^OTT 3 



eXP( -2^ 



(C-4) 



(C-5) 



we obtain localized a-motion such as 



gaussian 



f a \ 1/4 - 
x Vo (/3+,a)Vo (fa, a 



exp 



a 9 
2 a 



+ (-1 



i-K/2 



exp 



a 



(C-6) 



where the normalization constant is given for the case with no substantial overlapping 
between two terms in the square bracket on the r.h.s. of Eq. (|C-6j) . Note that due to 
the application of the selection rule K db v = Am on cos vol, the resultant a-function 
(the second line of Eq. (|C-6p ) is periodic with period tt, and especially for K = An 
the function is with period 7r/2. Our expression for the variable a in Eq. (|CJ -6|) is 
given for the region — tt/A < a < Sir /A, so that the next gaussian peak at a = tt can 
be omitted. 

The normal modes for the /?+- and /3_-motions are named as butterfly and anti- 
butterfly modes, respectively, the quantum numbers of which are (n+, n_). Due to the 
selection rule (— l) n ++ n - = (_i )^ ; the lowest butterfly and anti-butterfly states with 
K = appear with (n+,n_) = (2,0) and (0,2), respectively. The physical butterfly 
motion corresponds to the configuration displayed in Fig. 2, where the motion of the 
axis z'2 is confined around downside with ai ~ tt and the vibrational motions with 
A fa ~ A fii, for example, with the quanta (2,0). When the configuration described 
with variables a, and Afa is transformed by I2 ■ («2, fa) — K a 2 + 7r > 7r ~~ $2) i n 
Eq. (|B-3p . (/3 + ,/3_) is transformed into (/?_,/?+), which means that the physical 
butterfly motion is also described with configurations with the motion of the axis z% 
confined around upside with a<i ~ and the vibrational motions with Afa ~ —A fa, 
with the quanta (0, 2) as the example. Those symmetric terms for (n+, n_) exchange 
are described in Eq. (|B-9p and Eq. (|C-2p , where wave functions for the lowest butterfly 
state hK2o{ot, fa) consist with one term with (2, 0) of a ~ ir/2 and another term 
with (0, 2) of a ~ 0; for example, for butterfly, 



(— 



1/4 



{-if/ 2 



exp 



+ ex P ( ~2 a ' 



with a zero-point oscillation assumed for the a-degree, while for anti-butterfly, 
/ a \l/4 r f at 7r \ 2 



(C-7) 



a 

V47T 



(_l)*/2 



exp 



a 



+ exp 



-a 



(C-f 
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We have solved the Schrodinger equation for the a-degree in §3, where we worked 
with the reduced potential for the butterfly quanta (2, 0) and obtained the solutions 
4>k(o()- A butterfly state appears with the lowest energy with a localization around 
a = it/2 as is seen in Fig. 7. The motion for the a-degree of an anti-butterfly state 
which corresponds to the quanta (2, 0), i.e., 4>k{o) for the second term of Eq. (jCJ-8|) is 
obtained with an excitation. In order to calculate the partial widths and the angular 
correlations for the butterfly states, we use a simple but a typical expression, in 
which 4>k{.oi) are described with dominant two coefficients Cq and C 2 in Eq. ()C-2p . 
Furthermore to obtain a typical expression, we impose a symmetry between the a- 
motions in the butterfly and anti-butterfly states. Then we have the expression for 
the butterfly state with K = even as 

h K20 (a,P+,P-) = -^[(-l) x / 2 (l - V2cos2aV 2 (/3+Vo(/3-) 

+(1 + V2cos2a) ¥ >o09+)^03-)], (C-9) 

where the oscillator <p , (j3± , a) are also simplified to be independent upon a with 
an averaged oscillator energy. The corresponding pair of the anti-butterfly state is 
given by 

h K02 (a, /?+,/?_) = [ {-l) K ' 2 {\ - V2cos2a)99 (/3 + )^(/3-) 

+ (1 + v / 2cos2a)<^ 2 (/3 + )^ (/3-)]. (C-10) 

As the dynamical solutions for the a-motions in the reduced potential with the 
quanta (2,0), <Pk(®) are simplified to be (1 — \/2cos2q) for the butterfly mode, 
while they are (1 + \/2cos2q) for the anti-butterfly mode. Note that the results 
of our dynamical calculations with the reduced potential are rather close to those 
typical expressions. 

Wobbling motion (K -mixed states) 

Following the discussion on K- mixed states in §4.1, we define a wave function 
for the wobbling motion. The mixing weights are given by a gaussian function of K, 
by which we superpose D J MK {6i) such as 



T.^fW D MK(e i ) = {^^di IK {e 2 )Y^c K e ^ (on) 



K K 

with 

~exp[-l(f 



C'k- { \/V^ f L 2\,t; j for K ~ even (c-i2) 
for K = odd, 



where \AK\ = 2 due to the nature of couplings due to the axial asymmetry. The term 
with K = is important for including the components of the elastic scattering. The 
coherent summation over Ck^ %K6a again gives us a gaussian function ^(^3) of 
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period ir due to K = even. We obtain, for < #3 < 2ir, 

w(e 3 ) -- 



2^ 



exp { - b -el) + exp { - b -(9 3 -rr) 2 } + exp { - -(9 3 - 2vr) 2 } 



(C-13) 

The wobbling basis function of Eq. (|C-lip would be applied with the internal 
modes xk(R, a, /3_) of Eq. (|B-5p . which receives the transformations 2j given 
in Eq. (|B3p . and is symmetrized as in Eq. (|B-8p . When we sum over D J MK {0i) of 
different K- values, relative phases of those functions should be chosen properly. The 
coefficients for the wobbling motion in Eq. ()C-lip is given for the axially asymmetric 
configurations, where the molecular axis with the largest moment of inertia is x' , 
i.e., I x i > I y ', such as for an equator-equator one illustrated in Fig. 8. However 
the butterfly configuration such as illustrated in Fig. 2 has I x > < I y i, because, with 
a\ = and a 2 = 7r, i.e., #3 = tt/2 and a = —ir/2 by definition, the x'-axis moves 
to the direction of the initial y'-axis. To recover the condition I x i > I y r, we need to 
reset the x'-axis on to the intial direction, which brings additional phases e~ ln ^ 2 = 
(—i) K on D J MK (Qi). The expressions in Eqs. (|C-6p ^ (|C-10p satisfy this relative phase 
convention. Of course, we are able to adopt the rotational equation of motion with 
I x i < I y i in §4.1, to obtain the alternative phase (— 1) K I 2 for the K-bases. 



Appendix D 

Moments of Inertia of Dinuclear Systems 



The expression of the inertia tensor of dinuclear systems is already given in 
Eq. ()2-8p . which is defined by the configuration referring to the molecular axes. We 
again write it here, for convenience, i.e., 

/ s = Ip{R) + *#(ai0i7i) hR'iaxPm) + *#(a 2 /3 2 72 ) J 2 #(a 2 /? 272 ), (D-l) 

where the first term of the r.h.s. denotes the moments of inertia of two-ion centers, 
and the second and third terms are individual contributions from the constituent nu- 
clei with rotation matrices R' (a.ifii'yi). The diagonal components In and I 22 of the 
inertia tensor i^(-R) are fJ-R 2 , the others being zero. The inertia tensors of the two 
constituent nuclei, I\ and i" 2 are defined in the coordinate frames of their principal 
axes. Then, they are diagonal, elements of which are determined by the excitation 
energies of the members of the ground rotational bands of the constituent nuclei. 
Except for the relative vector of the two- ion centers, the whole dinuclear configura- 
tion is determined by the orientations of the principal axes of the constituent nuclei, 
due to Euler rotations fi[{ai, fy, 7j)^m(^1j #2)- In them, the first rotation Qm^i,^) 
is concerned about the molecular axes, and the second rotation f2' i (ai/3i r Yi) is that 
of each constituent nucleus referring to the molecular axes. Thus the moments of 
inertia about the molecular axes are obtained with the rotation matrices R '(o!iA7i)- 
We estimate magnitudes of the moments of inertia from the shape of the molec- 
ular configuration displayed in Fig. 8. Inserting a\ = a 2 = and f3\ = /3 2 = tt/2, we 
obtain the diagonal elements of I s to be 
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I x = fiR 2 + I a + I b , 

I y = iiR 2 + I a + Ib, (D-2) 
I z = I a + Ib, 

with the nondiagonal elements being zero. I a etc. denote the moments of inertia of 
the constituent nuclei, individually in their principal axes, i.e., the diagonal elements 
(Ix,l2,l3) of I\ are written as {I a, I a, la), and those of I2 as (Ib, Ib, h), and their 
nondiagonal elements are zero. Note that for the states of the ground rotational 
band of the 28 Si nucleus, due to the axial symmetry, I\ = I2 is assumed and further 
I a = h = is adopted in §2. 

The value of the moment of inertia I for the 28 Si ground band, is determined 
from the excitation energy E x = 1.78MeV of the 2f state of the 28 Si nucleus, i.e., 
by the relation, 

g/(/ + l) = £ x , (D-3) 

/ being the spin value (1 = 2), and the value of lis used in the numerical calculations 
in §3. Note that I A = Ib = Z is denoted by / in Eq. flS-H}- 

On the other hand, moments of inertia of rigid bodies have been often investi- 
gated in the study of rotational spectraP^ In §4.1, with respect to the nuclear shape 
of the whole system, we adopt I a and Ib of nonzero value, where induced deformation 
is expected. Moments of inertia can be obtained by integrating over nuclear volume, 
such as 

h = f p(r)(r 2 - x 2 )dV, (D-4) 

where Xj denote the coordinates in the principal axes, and p(r) is a nuclear density 
distribution. The density profile of 28 Si appears in the calculations of the folding 
potential with DDM3Y force in §2.3, and its parameters are given there. With 
induced deformation, I\ is not necessary to be equal to I2, but we take the value 
of I\ = I2 tentatively, due to the axial symmetry of the density profile. Due to the 
oblate shape of the density profile, we obtain the vales of 1% = I2 < I3, the values 
of which are (164, 164, 222) in the unit of M n fm 2 , with M n being the nucleon mass. 
Compared with the value of moment of inertia, 70M n fm 2 estimated by Eq. ()D-3j) . 
the value 164M n fm 2 is about two times larger than that, which is well known for 
rotational spectra of nuclei.^ Since the moments of inertia of rigid body are too 
large, we renormalize the values of moments of inertia obtained by Eq. (|D-4p . to 
be consistent with the excitation energy of the 2f state of the 28 Si nucleus. This 
means that we multiply a factor 0.42 on the moments of inertia of Eq. (|D-4j) . Thus a 
value I a = Ib = 93M n fm 2 is adopted, and Eqs. ()D-2p give the values for moments of 
inertia of the whole system. Note that I\ = I2 is broken with induced deformation 
of the 28 Si nucleus, but we have no information about those changes of the moments 
of inertia. Hence we adopt the same value 70M n fm 2 for Ia with the assumption 
I\ = I2, as well as for Ib- 
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As for dinuclear configurations in §4.2, since one of the constituent nuclei is 
spherical, i.e., 1 2 = in Eq. (|D-ip . we obtain the expressions of moments of inertia 
by simply inserting Ib = lb = into Eq. (|D-2p . Note that those expressions are 
tailored for the configurations displayed in Fig. 11(b) and Fig. 12, where Ia and I a 
are denoted as Jo and I3, respectively, because of no 1 2- For case 1, I a = is adopted 
with the axially-symmetric deformation, and for case 2, I a ^ is adopted with the 
axially-asymmetric one, in which the ratio I a /lA can be estimated by Eq. (|D-4|) . 
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